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ABSTRACT 



This report describes and evaluates a practical computer based 
method for translating data concerning: 

1. the location of each school to be serviced by a bus fleet, 

2. the locations and numbers of students to be transported to 
each school, 

3. the time interval during Which the students are to be 
transported, and 

4. the available bus facilities 

into a set of bus routes which specify school- to- school sequencing of 
each bus and the stop-to-stop route to be followed in traveling to 
every school. Each route is designed in such a way that the bus 
capacity and student riding time constraints are satisfied while 
attempting not only to minimize the total bus travel time (including 
running empty) for a school but also to minimize the number of routes 
required to service all the stops associated with the school. The 
mathematical models developed were programmed in FORTRAN IV for use on 
a CDC 6400 computer and were applied to four schools in the Williamsville 
New York Central School District. An efficient routing system 
involving 6 possible bus route origins and 96 stops was developed for one 
of these schools in 61 seconds on a CDC 6400 computer. 
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I. INTRODUCTION 



School districts are aware of the power of the digital computer as 
a tool for both reducing cost and improving service in the management 
of their educational programs. The computer has already proven its worth 
in financial accounting, in personnel administration, in class scheduling, 
and in planning school construction. 

One area of school administration which has not yet been adequately 
served by the capability of the computer is the management of the 
transportation system. Since the cost of procurement, maintenance, and 
operation of a bus fleet requires a large portion of a school district's 
budget, the director of transportation is expected to minimize these 
costs while simultaneously providing an acceptable level of service. 

At present, most schools prepare bus routes and schedules manually 
by using a large map of the district and a listing of the school census. 
Since a single school may have as many as two hundred bus stops and a 
district as many as twenty-five schools, this procedure is not only time 
consuming but requires an excessive amount of administrative talent 
which could be better utilized in other endeavors. Moreover, the quality 
of the bus routes prepared by hand is a function of the scheduler's 
experience, i.e. the best routes are usually prepared by the most ex- 
perienced schedulers. 

In addition to the bus scheduling problem, school administrators 
are also interested in being able to evaluate the sensitivity of current 
fleet operations to various changes. For example, it would be useful to 
be able to easily examine how population fluctuations, a change in bus 
fleet size, a change in school boundaries, or a change in policy con- 
cerning the maximum allowable walking distance would affect bus operations. 
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Therefore, a great need exists for some means for generating school bus 
routes and schedules efficiently. 

This dissertation describes and evaluates a computer based 
methodology for translating data concerning: 

1. the identification and location of each school which is to 
be serviced by the bus fleet, 

2. the location and number of students to be transported 
to each school, 

3. the time interval during which students are to be 
'transported, and 

4. the available bus facilities 

into a set of bus routes which specify school-to-school sequencing of 
each bus and the stop-to-stop route to be followed in traveling to the 
school. Each route is designed in such a manner that the bus capacity 
and maximum allowable student riding time constraints are satisfied 
while attempting to not only minimize the total bus travel time (includ- 
ing running empty) for a school but also to minimize the number of routes 
required to service all the stops associated with the school. 

The output of the computer program based upon this methodology is 
the detailed information needed to prepare specific bus schedules and 
bus passes for each student. In addition it provides a means for de- 
termining overall measures of schedule performance as total travel time, 
average riding time, average bus load, etc. These performance measures 
coupled with the inexpensive and rapid computer development of a complete 
set of bus routes for a school district make feasible a quantitative 
evaluation of the effect of any administrative policy changes on the 
transportation system. 



The computer model makes operationally feasible the generation of 
bus routes and schedules immediately prior to the start of the school 
year when knowledge of transportation demands is most accurate. More- 
over, it provides a tool for evaluating alternative new school sites in 
a district experiencing significant population growth. 



II. LITERATURE REVIEW 
A. School Bus Routing Problem 

A review of the literature concerned with school transportation 
reveals that most of the studies conducted in this area deal with 
standards for school buses, methods of allocating funds to school dis- 
tricts for transportation purposes, surveys of the status and philosophy 
of school transportation in various geographic regions, and examination 
of legislation affecting school transportation. The absence of extensive 
literature dealing with the school bus routing problem indicates that 
either interest in this problem area is relatively recent or that a large 
portion of the attempts at developing a general method for school bus 

routing have proved to be fruitless and have therefore not been reported. 

12 

Dantzig, Fulkerson, and Johnson gave the allusion that the problem 
of using a digital computer to design school bus routes was simple and 
straightforward. Intrigued by this, Boyer"* attempted to solve a five 
stop school bus problem by applying the simplex algorithm. Because the 
application of a linear programming technique to the problem seemed to 
be highly impractical, he developed the Sequential Steps Method^. This 
procedure is based upon the premise that all students should be tran- 
sported the shortest possible distance to school. Each route starts at 
the bus stop which is the farthest distance from the school and proceeds 
along the shortest path to the school in such a manner that no isolated 

bus stops are created and that a bus stop is assigned to only one route. 

34 

Thompson tested the effectiveness of the Sequential Steps Method 
in determining school bus routes for a hypothetical school. He learned 
that although student riding time was lowered, this method tended to 
generate more routes for a school than other manual methods favored by 
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the transportation personnel participating in the study. Although this 
method was sequential, it was not amenable to computer programming. 
Moreover, it made no provision for imposing bus capacity and passenger 
riding time constraints. 

Boyer^’^ described a procedure for designing school bus routes in 
which a set of bus routes or a route system is developed manually. A 
computer is then used to list all possible ways in which the stops of a 
bus route can be visited. By inspection, tne best permutation for each 
individual route of the route system is determined. After several 
arbitrary route systems have been analyzed in this manner, the best 
route system is selected. 

Since this method requires a great deal of manual work and personal 
judgement, it offers little improvement over the current manual methods 
of designing bus routes and schedules. Moreover, this method makes no 

provision for imposing a passenger riding time constraint. 

35 

Tillman applied the technique of dynamic programming to the school 
bus routing problem in which all routes start and end at the school. In 
this ^formulation stage j was the number of the bus being loaded or the 
number of the route being designed, the decision variable was the number 
of stops made by bus j or the number of stops assigned to route j, the 
stage input was the number of stops not yet serviced at stage j, and the 
return obtained at each stage was the minimum number of miles traveled 
by bus j. The objective was to minimize the sum of the returns subject 
to a bus capacity constraint. Using this method, an optimum solution 
was obtained for a problem involving five bus stops, three buses, and 
forty students. 

Although dynamic programming guarantees an optimal solution and 



allows the interstop travel time matrix to be asymmetric, it is an 
impractical method because of the extremely large number of calculations 
which have to be performed for even relatively small school bus routing 
problems. 

Newton and Thomas^®’ described a practical method for generating 
school bus routes and schedules by computer. Given the matrix of inter- 
stop travel times, which may be asymmetric, bus routing is accomplished 
by a two step procedure. First, a single near-optimal route which starts 
at the school, visits every stop once, and terminates at the school is 
determined. This route, the solution of the traveling-salesman problem 
associated with the given set of bus stops, is then partitioned into 
individual bus routes which satisfy bus capacity, bus loading policy, 
and passenger riding time constraints. The order of the route determined 
in step one is preserved during the partitioning process and all routes 
originate and terminate at the school. This heuristic procedure has 
been used to solve an eighty stop problem in approximately six minutes 
on a 7090 computer. 

Davis ^ described a branch and bound algorithm for solving the 
school bus routing problem in which the interstop travel time matrix is 
symmetric, bus capacity and passenger riding time constraints are imposed, 
and all routes start and end at the school. This procedure partitioned 
the set of all possible routes into mutually exclusive subsets of routes 
by solving a series of transportation problems. The cost of each sub- 
set was the total traveling time required by the optimum solution to the 
transportation problem associated with the subset. 

Although this method is amenable to computer programming and guaran- 
tees an optimal solution, it proved to be an impractical one. For 
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example, an attempt to solve a thirteen stop bus problem, whose interstop 

travel time matrix was symmetric, subject only to a bus capacity con- 
14 

straint from was abandoned when no feasible solution was obtained after 
approximately eighteen minutes CDC 6500 time had been spent on the 
problem. 

B. Delivery Problem 

Considerably more information is available on the closely related 
delivery problem*■*®*®*'^ , The delivery problem is concerned 

with the determination of routes for a vehicle, initially located at a 
depot, which visits a number of delivery or pickup points and returns to 
the depot. Since the capacity of the vehicle is less than the total 
quantity of goods which must be transported, several trips must be made. 

The delivery problem usually is not subject to a traveling time or dis- 
tance constraint whereas the school bus routing problem is nearly always 
constrained by both bus capacity and maximum allowable student riding 
time. 

Balinski and Quandt* formulated the delivery problem as an integer 

programming problem. Although this method guarantees an Lptimal solution, 

it is undesirable because of the large number of variables and constraints 

required to express a delivery problem involving relatively few stops. 

Moreover, the available integer programming algorithms are often unable 

to achieve solutions to even moderately large problems even though 

theoretically they should always determine the optimum solution. 

21 

Hayes applied the branch and bound method of Little, Murty, Sweeney, 
27 

and Karel to the delivery problem. This procedure partitioned the set 
of all possible routes into mutually exclusive subsets by either assign- 
ing or hot assigning a particular link to a route. The cost of each 
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subset was the minimum possible total length of any set of routes contain 
ing the link assignment associated with the subset. 

Although this method guarantees an optimal solution and is amenable 

to computer programming, it also is an impractical method because a large 

number of time consuming operations must be performed. For example, an 

14 

attempt to solve the thirteen stop problem from was abandoned when no 
solution was obtained after one hour was spent on the Control Data G-21 
computer. 

14 

Dantzig and Ramser described a heuristic procedure for solving the 
delivery problem when all trucks have the same capacity. After the 
delivery points are arranged in numerically ascending order with respect 
to demand and the number of stages required to achieve solution has been 
calculated, this algorithm synthesizes routes by a stage-wise aggregation 
of the delivery points, i.e. in stage one, pairs of points are joined, in 
stage K, groups of K points are joined to other groups of K points. At 
each stage the points are combined in such a manner that truck capacity 
is not exceeded and the sum of the interstop distances is minimized. 

. Although this method is sequential and. amenable to computer pro- 
gramming, personal judgement may be required in adjusting the final 
solution which may not be uniquely defined. Moreover, this procedure 
places more emphasis on insuring that the trucks are loaded to capacity 
than on minimizing the total distance traveled. It also appears that 
this algorithm would require considerable alteration to handle the 
asymmetric case because direction is not considered when the interstop 
distance matrix is searched to determine the minimum entry. This method 
would probably be quite time consuming in solving problems which involve 
more than 20-30 stops. Each entry of the distance table in stage K, 
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where K ■ 2,3,...,N and N is the number of stages required, contains the 

shortest route which starts at the depot, visits 2 points and returns to 

the depot. Therefore, many traveling-salesman problems must be solved at 

each stage. For example, if after 2 stages 100 stops are aggregated into 

25 groups of 4 stops each, then the distance table for stage 3 requires 

the solution of 300 traveling-salesman problems of 9 stops. 

o 

Clarke and Wright described a heuristic procedure for solving the 
delivery problem when all vehicles do not have the same capacity. After 
the demand points are arranged in numerically ascending order with 
respect to distance from the origin, the algorithm assigns each stop to 
one vehicle or route. Then each pair of points is examined to determine 
the savings which would result if they were linked together on the same 
route instead of being assigned to different routes. A savings calcu- 
lation is made only for those points, currently linked to the origin, 
which could be assigned to the same route without violating vehicle 
capacity constraints. The two points associated with the maximum savings 
are then assigned to the same route. After repeating the procedure 
until no further savings can be calculated, the resulting set of routes 
is the solution to the problem. 

This method is amenable to computer programming and appears to be 
efficient for solving the delivery problems in the literature. However, 
it cannot be readily extended to handle the case of the asymmetric 
inters top distance matrix because direction is not considered in calcu- 
lating the savings or the criterion for accepting or rejecting the 
linkage of two points on the same route. Moreover, the largest problem 

for which solution success with this method has been reported is one 

18 

involving thirty- two stops . 
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Cochran^’ ^ made two modifications to the algorithm of Clarke and 

s 

Wright . The first change permits the reassignment of vehicles to in- 
dividual routes each time a vehicle becomes available as a result of the 
linkage of two points on the same route. The second change allows an 
upper bound to be placed upon the length of any route. Although the 
first modification insures that the vehicles will be more fully utilized 
and the second modification is a useful one, neither the probability of 
reaching optimality nor the efficiency of the algorithm is increased by 

their inclusion. 

8 

Braun described a simulation approach to the delivery problem. 
First, delivery points are randomly assigned to an individual route sub- 
ject to the vehicle capacity constraint. Each route is then improved 
by applying a traveling-salesman algorithm. After a specified number of 
sets of routes are developed in this manner, the set which covers the 
least number of miles is selected as the solution. 

Although this procedure is simple and amenable to computer pro- 
gramming, it proved to be less efficient than other available methods 
with respect to the quality of the routes produced and the amount of 
computer time required to achieve solution. Moreover, the quality of 
the sets of routes tends to become poorer as the number of delivery 
points increases. 

After abandoning the branch and bound method for solving the 
21 

delivery problem, Hayes developed a computerized version of a procedure 
which an experienced dispatcher might employ for routing trucks. Since 
this method assumes that the warehouse or depot is located near the 
center of the scatter of customers and performs poorly when this assump- 
tion is not satisfied, it does not appear to be superior to the available 



heuristic methods. 

18 

Gaskell reported on his experiments with the method of Clarke and 
9 

Wright . He devised five functions for calculating the savings which 
would result from linking two points on the same route. After applying 
each of them to a set of delivery problems, he concluded that none was 
uniformly better and that the function used by Clarke and Wright to 
calculate savings was a reasonable one. 

Although more work has been done on the delivery problem, success- 
ful solution has been reported only for relatively small problems with 
no indication of apparent near future breakthroughs for larger problems 
involving asymmetric interstop time/distance matrices, routes whose 
origin and terminus do not coincide, and restrictions upon the length 
of a route, vehicle capacity, and the number of routes used to service 
all the 8 tops. Indeed, solution success on the school bus routing pro- 
blem may provide a means of solving the delivery problem. 
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III. PROBLEM FORMULATION 



A. Statement 

Thus far, the school bus routing problem has been solved for an 
environment which has been greatly simplified by the imposition of 
various assumptions. Most solutions not only exclude constraints on 
some of the route characteristics but also assume that the origin and 
terminus of a bus route coincide, whereas in reality they are not 
necessarily the same. In fact, there often may be several possible 
origins for the set of bus routes servicing a particular school. There- 
fore, it would be desirable to be able to judiciously select the proper 
combination of origins from the set of possible origins for the system 
of bus routes servicing a particular school. 

With these thoughts in mind, this dissertation is concerned with 
the determination of all bus routes for a school district. Buses are 
routed from school to school picking up students as they travel. A 
heuristic procedure has been developed and programmed in FORTRAN IV to 
generate efficient bus routes when the school bus routing problem has 
the following definition: 

Given: 
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1. the number of time periods used by the school district for 
bussing 

2. for each time period 

a. the number and location of schools to be serviced 

b. the number of possible origins for the bus routes 

c. the identification of each origin 

d. the number of buses available at each origin at the 
beginning of the period 

12 
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3. for each school to be serviced during the same time period 

a. the identification of the school, 

b. the identification of each stop to be visited, 

c. the number of students assigned to each stop, 

d. the matrix pf Interstop travel times between each 
possible origin and every stop assigned to the 
school, 

e. the matrix of interstop travel times for all pairs 
of stops except those links involving an origin, 

f. the bus capacity, 

g. the maximum allowable student riding time, and 

h. the criterion for accepting a set of feasible routes 
as the quasi- optimal solution. 

Determine for each school: 

1. the set of bus routes and schedules required to provide 
transportation for all students either to school or from 
school and 

2. a lower bound on the total number of time units required 
to traverse all the bus routes 

such that: 

1. no more than the absolute minimum number of routes required 
to transport the students plus one route will be used for 
any school, 

2. the bus capacity and maximum allowable student riding time 
constraints will be satisfied, 

. the criterion for accepting a set of feasible routes as 
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the quasi-optimal solution will be satisfied, and 

4. optiraizaticfn will be with respect to minimizing the total 
traveling time for a set of routes. 

It is assumed that; 

1. all buses assigned to a particular school have the same 
capacity, 

2. all routes for a particular school will be subject to the 
same maximum allowable student riding time constraint, 

3. the matrix of interstop travel times may be asymmetric, 

4. a stop is assigned to only one route, 

5. the traveling time from any bus stop to the terminus Is 
less than the maximum allowable student riding time, 

6. the origin and terminus of any bus route do not necessarily 
coincide, 

7. a bus services only one route for a particular school, and 

8. the number of buses available at the beginning of period i 
is adequate to meet the needs of all students to be serviced 
during period i, i.e. it is possible to service all the 
routes simultaneously. 

B. Discussion of the Problem Statement 

Most districts set the opening and closing times of the schools 
serviced by the same bus fleet so that all students can be transported 
during two or three non- over lapping bussing periods whose total elapsed 
time is approximately two and one-half hours. Schools maintaining 
common hours of operation are serviced during the same bussing period. 

The data required to design the routes for each school can be 



developed either manually or by computer from the map of the area and 
. the associated census. tract. Each stop serviced by the school Is 
usually Identified by both name and number In order to facilitate com- 
munication between the administrators and the bus drivers. The elements 
of the Interstop travel time matrix are calculated from a map of the 
area so that they include the effects of both distance and expected driving 
conditions in traveling between every pair of points. Often in an effort 
to reduce the size of the interstop travel time matrix, bus stops ad- 
jacent to each other on the same side of the road are combined into a 
single stop. The information required to assign students to a bus stop 
can be extracted from the census tract. Usually upper limits are im- 
posed upon the number of students assigned to a bus stop in order to 
reduce the noise level and the possibility of landscape damage in 
residential areas. 

Optimization is with respect to not only minimizing the total time 
required to traverse a set of routes but also with respect to minimizing 
the number of routes required to transport all the students of a 
particular school. Both of these factors contribute heavily toward the 
operational costs of maintaining a bus fleet. Minimizing the number of 
routes insures that the buses are being utilized to full capacity and 
also tends to reduce the total traveling time associated with a set of 
routes by eliminating some of the links between an origin and the first 
stop of any route and the links between the last stop of any route and 
the terminus. Since schools usually place more emphasis on minimizing 
the number of routes required to transport all the students than on 
minimizing the total traveling time for a set of: routes, the maximum 
allowable number of routes for a school will be considered to be a 
constraint. 

15 
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The minimum number of routes required to service all the stops 
assigned to a school is the smallest integer which is greater than or 
equal to the quotient of the total number of students to be transported 
and the bus capacity. Since a riding time constraint is imposed on each 
route and since all the students assigned to a particular bus stop have 
to be picked up by the same bus in order to avoid confusion, it may be 
impossible to service all the bus stops by the minimum number of routes. 
Therefore, a set of routes is considered feasible if each route satisfies 
the bus capacity and student riding time constraints and if the set con- 
tains at most one more than the minimum number of routes required to 
transport all the students. Thus, a feasible set of routes contains 
either the minimum number of routes or one more than the minimum number 
of routes. 

Because the set of bus routes for any school is designed by a 
heuristic procedure, some method must be devised to determine when a set 
of routes, acceptable to the school, has been developed. One possible 
criterion is to accept the best set of feasible routes available after 
the procedure has been executed a specified number of times. A second 
possible criterion is to accept the first set of feasible routes deve- 
loped whose total travel time is less than the product of a given factor, 
greater than one, and the lower bound upon the total time required to 
traverse all the routes for a school. Another alternative is to accept 
as the quasi-optimal solution the best set of feasible routes available 
at the time at which either of the criteria is first satisfied. 

The lower bound upon the total time required to traverse a set of 
routes for a school is the minimum length of time in which all the bus 
stops could be serviced by the number of routes included in the set. 
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In calculating the lower bound, it is assumed that all connections bet- 
ween' the origin and the first stop of a route, all the connections 
between the bus stops, and all the connections between the last stop of 
a route end the terminus are made in the least time consuming manner. 
Since the lower bound does not include the effects of the bus capacity 
and passenger riding time constraints, the probability of developing a 
set of feasible routes whose total traveling time equals the lower bound 
is very low. However, even though the lower bound may be unattainable, 
it provides some measure for assessing the quality of the set of routes 
developed. 

An examination of the set of assumptions under which the problem is 
to be solved reveals that the environment has not been oversimplified. 
The bus capacity is based upon the expected girth of children at various 
age levels, the size of the bus used, and the degree of bus utilization 
required by the school. The maximum student riding time is dependent 
upon the size of the area, the scatter of the bus stops and the local or 
state laws regulating student riding time. Usually, districts consider 
the bus capacity and the maximum riding time constraints to be fixed for 
a particular school. 

The interstop travel time matrix is characteristically asymmetric 
for reasons such as; restrictions imposed upon the crossing of busy 
streets by children, one-way streets, limited access highways, and re- 
strictions imposed on vehicle turns at intersections. Moreover, most 
schools maintain a bus loading policy which requires that all students 
assigned to a stop by picked up by the same bus in order to avoid con- 
fusion. This restriction implies that a stop is assigned to only one 
route . 
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The traveling time from any bus stop to the terminus has to be less 
than the maximum allowable, student riding time in order to insure that 
no bus stop b- isolated. The assumptions that a bus services only one 
route for a school and that it is possible to service all the routes 
simultaneously are imposed in order to facilitate the computational pro- 
cedure used to develop all the bus routes for a district and do not 
simplify the environment. Thus, this formulation of the school bus 
routing problem is considered to be a realistic and a reasonable one. 
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C. Mathematical Model 



The bus routing problem for any school, as previously defined, can 
be expressed as a zero-one integer programming problem. 



j ■ terminus of any link 
k « number of any possible bus origin 
m «* route number 

K “ number of possible different bus origins for a set of 
routes 

M = number of routes allowed for the school 
N = total number of bus stops (includes the origin and 
the terminus) assigned to the school 



stop 1 (k) = bus origin number k 

stop N = terminus of any route or the school 
t(i,j) - number of time units required to travel between 
bus stop i and bus stop j 
L(j) = number of students assigned to bus stop j 
B(k) = number of buses available at origin k 
C - bus capacity 

R = maximum allowable student riding time 
x(i,i-m^ a C 1, if link (i,j) is assigned to route m 



The subscript denoting the particular school has been omitted in order 
to simplify the nomenclature. The words bus and route are used inter- 
changeably. 

The problem can be stated as follows: 

Find variables x(i,j,m) for all combinations of i, j, k and m where 



Let 



i «* origin of any link 




0, otherwise 



i = l(k),2,..., (N-l) 
j = 2,3,... ,N 

k = 1,2 K 

m 1,2, . . . ,M 



such that the objective function 
N K M 

X Z Z Ct(l(k),j)][x(l(k),j,m)] 

j=2 k=l m=l 



+ 



(N-l) 

i 



i=2 



N M 

£ X Ct(i,J)][x(l,J f n)] 
j=2 m=l 



is minimized while restrictions one through eight, stated beiow, are 
satisfied. 

Optimization, as previously stated, is with respect to not only 
minimizing the total time required to traverse a set of routes but also 
with respect to minimizing the number of routes required to transport 
all the students of the school. Since it is difficult to work with 
two objective functions, and since most schools place greater emphasis 
on the minimization of the number of routes than on the minimization 
of the total traveling time, the number of routes required to transport 
all the students of a school will be considered a constraint and the 
objective function will involve the total traveling time only. The 
first term of the objective function is the total time spent in 
traveling from the origin to the first stop of any route; the second 
term is the total time spent in traveling between the other pairs of 
points. 




20 



The first constraint, 

I 

KM / (N-l) M 

—l -“I — 1 •— I 

L L + 2, 2 J “ 1 " for i = 2,3,,.., (N-l), 

k“l m=l i=2 m=l 

insures that any bus stop j, where j is not the origin or terminus of 
any route, will be the terminus of exactly one link on one route. There 
will be (N-2) constraints of this type. 

The second constraint, 

N M 

£ 2, ■ 1 > for 1 ■ 2,3,..., (N-l)., 

j=2 m=l 

insures that bus stop i, where i is not the origin or terminus of any 
route, will be the origin of exactly one link on one route. There will 
be (N-2) constraints of this type. 

The third constraint, 

(N-l) K M 

E Z 2j x ( 1 ( k )»j»'ft) = M, 

j=2 k=l m=l 

insures that only M routes are used by counting the number of links 
between the origin of each route and the first stop serviced by the 
route. 

The fourth constraint, 

(N-l) M 

X l x(i,N,m) = M, 
i=2 m=l 

also insures that only M routes are used by counting the number of 

O 
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links between the last stop serviced by a route and the terminus or 
school. 

The fifth constraint, 

(IM.) M 

2, *(l(k),j,m) £ B(k), for k « 1,2, .. . ,K, 

j=2 m=l 

insures that the number of routes starting at a particular origin does 
not exceed the number of buses available at that origin for the school. 
There will be K constraints of this type. 

The sixth constraint, 

(N-l) K (N-l) (N-l) 

2, Lx(l(k), j,m)][L(j)] + 2, Z i- x (i» * C, 

j=2 k=l i=2 j=2 

for m = 1,2,...,M, insures that the bus capacity will not be exceeded 
by any route. The first term of the constraint counts the number of 
students picked up at the first bus stop serviced by the route and the 
second term counts the students picked up at the other stops assigned 
to the route. There will be M constraints of this type. 

The seventh constraint, 

(N-l) N 

— i —l 

2, l C*<i»J>3C*<i» * R. for m = 1,2, ...,M, 
i=2 j=2 

insures that no route will exceed the max imum allowable student riding 
time where the riding time is counted from the first pick-up point. 
There will be M constraints of this type. 

The eighth constraint, 

(r-1) 

■*» 

2 , x ( i p» i p +i» in ) + * (r-l), for m = 1,2, ...,M, 

P=1 
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where i^,i 2 *r ran S es over a H permutations of subsets of the bus 

stops {2,3, . . . , (N-l)J of size r, and 2 £r £ (N-l), prevents the for- 
mation of a loop, a route which starts and ends at the same point. 

(N-2) z 

There will be 51 / P^ N constraints of this type where P denotes the 

^ k x 



k=2 



permutation of z things taken x at a time. 



Since 



(N-2) 

l 



k=2 



(N-2) 

P 

k 



is 



greater than 2(N-2).', there will be more than 2M(N-2).' constraints of 
this type. None of the permutations involve the origin or terminus 
because no link is allowed to start at the terminus and no link is 
allowed to end at the origin in the case of the general route whose 
origin and terminus do not coincide. This is accomplished by setting 
the travel times associated with these links equal to infinity, i.e. 
the elements of column one and row N of the inters top travel time 
matrix are assigned large values. To illustrate the manner in which 
the eighth constraint prevents the formation of a loop, let r = 2, 
m = 4, i^= 3, i 2 = i f - 5, x<3,5,4) = 1, x(5,3,4) = 1 and the rest of 

the variables, x(i,j,4) equal zero. In this example, route four starts 
at stop three, proceeds to stop five, and then returns to stop three. 
According to the eighth constraint, x(3,5,4) + x(5,3,4) = 2 = r and 
route four is a loop. 

This zero-one integer programming model would have to be solved 
twice. First, it would be solved for M equal to the minimum number of 
routes required to transport all the students of the school and then 
for M equal to one route plus the minimum number of routes required by 



O 

ERLC 



23 



the school. 



Since i may assume any of (K+N-2) values, j any of (N-l) values 
and m any of M values, there are (K+N-2) (N-l) (M) variables, x(i,j,m), 
in this problem. Moreover, there are more than 2M(N-2).' constraints. 

A normal bus routing problem for one school involving ten origins, 
twenty routes and one hundred bus stops would require 213840 variables 
and more than 40(98).' constraints. Since the available integer pro- 
gramming algorithms are clearly unable to handle problems of this 
magnitude, a heuristic procedure is the only recourse for developing 
all the bus routes for a school district. 
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IV. METHOD OF SOLUTION 



A. General Procedure 

Bus routing for a school district is accomplished by suboptimiza- 
tion from period to period. Bus routing for a period is accomplished 
by suboptimization over the sets of routes developed for each of the 
schools serviced during the period. Quasi-optimal bus routing for an 
individual school is accomplished by an algorithm based selective 
enumeration procedure. A detailed description of each step of this 
enumeration procedure is given in the following sections. It can be 
summarized as follows: 

1. Assuming that the maximum allowable number of routes will be 
used by all schools, determine for every school the number 
of bus routes which should start from each origin supplying 
buses for the period during which the school is being 
serviced so that the total estimated traveling time required 
by all the routes for the period is minimized. 

2. For each school, select the location which serves as the 
origin for the greatest number of routes as determined in 
step one. This point will be called the ''super-origin" for 
the school. 

3. Determine by either using the Nearest City Approach or 
Algorithm A, a trial route which starts at the super-origin, 
visits every stop once and terminates at the school. The 
method used to develop this trial route is dependent upon 
the number of times step three has been executed. 

4. Partition the single route determined in step three into 
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individual routes which satisfy bus capacity and passenger 
riding time constraints. The order of the stops determined 
in step three is preserved during the partitioning process. 

5. If the set of routes contains no more than the maximum 
allowable number of routes, then proceed to step six. 
Otherwise, return to step three to generate a different 
trial route. 

6. Improve each of the individual routes by Algorithm A, a 
modified traveling- salesman type algorithm. This step 
attempts to reduce the traveling time required by the 
individual route while preserving the assignment of stops 
to the route made by the partitioning procedure. 

7. Determine the current best set of routes developed by the 
procedure. 

8. If the set of routes obtained in step seven satisfies the 
acceptance criterion specified by the school, then proceed 
to step nine. Otherwise, repeat steps three through eight 
until either the acceptance criterion is satisfied or the 
algorithms used to generate the trial route in step three 
are exhausted, i.e. they are unable to generate another 
different trial route. 

9. If the results of step one specify that all the routes for 
the school should start at the super-origin , then this 
acceptable set of routes is considered to be the quasi- 
optimum solution. Otherwise proceed to step ten. 

10. Allocate the remaining origins specified by step one in 

such a way that the additional number of time unit 3 traveled 
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will be minimum. Then, this modified set of acceptable 
routes is considered to be the quasi-optimum set of routes 
for the school. 

After bus routing has been completed, the schedule or timetable 
for each individual route is calculated. Each schedule gives the bus 
load and the arrival time at each stop serviced by the route. The 
time required to load the bus at each stop and to unload the bus at 
the school is not included in calculating the timetable. In addition, 
a lower bound upon the total number of time units required to traverse 
all the routes required by the school is calculated. 

B. Selection of Bus Route Origins 

At the beginning of the first bussing period all buses are 
located at the garage(s) maintained by the school district. For all 
other bussing periods the available buses are initially located at the 
schools serviced during the previous period. It is assumed that the 
number of buses available at the beginning of any bussing period is 
adequate to meet the needs of all students to be transported during 
the period, i.e. it is possible to service all the routes for the 
period simultaneously. For each school j serviced during the same 
period, the determination of the number of bus routes which should 
start at each origin i can be expressed as a transportation problem. 
Let nl = number of origins supplying buses for the period 

n2 = number of schools serviced during the period 
x(i, j) - number of routes starting at origin i and ending 
at school j 

c(i,j) = average estimated time required to traverse any 



route starting at origin i and ending at school j 

A(i) = number of buses available at origin i at the 

beginning of the period 

R(j) m number of routes required by school j 

B(j) = bus capacity specified by school j 

The problem can be stated as follows: 

Find variables x(i,j) for all combinations of i and j, where 

i = l,2,...,nl and j = l,2,...,n2, such that the objective function 
nl n2 

J Y Cc(i,j)] [x(i,j)] 
i=l >1 

is minimized and restrictions one through four, stated below, are 

satisfied. Optimization is with respect to minimizing the total 

estimated traveling time required by the n2 sets of routes associated 

with the period. 

The first constraint, 
nl 

y x(i,j) = R(j), for j = 1,2,..., n2, 

i=l 

insures that all the routes required by school j are assigned to an 

origln c There will be n2 constraints of this type. 

The second constraint, 
n2 

£ x(i,j) = A(i) , for i = l,2,...,nl, 
j»l 

insures that that number of bus routes starting at origin i equals 
the number of buses available at origin i. It is assumed that a bus 
services only one route for a school. There will be nl constraints 
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of this type 



The third constraint, 
x(i,j) £ 0, an integer, for all i and j, 
is self-explanatory. 

The fourth constraint. 



nl 



n2 



V A(i) - £ R(j) 

j»l 



i=l 



insures that a feasible solution exists. It can be verified by 

n2 nl 



observation. Since the first constraint implies ) / x(i,j) 

fil 



n2 n2 nl 

y rq) and the second constraint implies Y ^ x(i,j) = 

>1 j=l i=l 

nl 

r» 

} A(i), then indeed the fourth constraint must be satisfied. 

i=l 

til n2 

Because it is assumed that / A(i) ^ y R(j) at the beginning 



i=l 



j=l 



of any bussing period, a fictitious school may have to be introduced 
to use the extra available buses in order to satisfy the fourth 
constraint. 

The parameters nl, n2, A(i) , and B(j) are given data for the 
problem. However, the parameters R(j) and c(i,j) for i = 1,2, ...,nl 
and j = l,2,...,n2 must be calculated. The R(j) are determined 



as follows: 

Let L(j,k) = number of students assigned to stop k 

associated with school j 
n(j) = number of stops assigned to school j 
(including the origin and terminus) 
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m(j) = minimum number of routes required to transport 
all students of school j 

The values of L(j,k) and n(j) where j = l,2,...,n2 and k = l,2,...,n(j) 
are given data for the problem. 

Tl(j) 

m(j) = r y LCjjk^/^BCj) J + 0.999999 truncated to the 
~ k=l 

nearest integer. 

R(j) = m(j) + 1, for j = 1,2 n2 

The parameter R(j) is used in the determination of the number of 
routes which should start at origin i for school j instead of m(j) 
because the bus capacity and student riding time constraints often 
make the generation of a set of feasible routes containing only 
m(j) routes unattainable. 

The c(i,j) are calculated as follows: 

Let LB(i,j) = lower bound on the total traveling time required by 
R(j) routes, all of which start at origin i and end 
at school j 

c(i j) = LB(i, j) / R(j) for i =1,2 nl and j = l,2,...n2 

The LB(i,j> are calculated under the assumption that all the 
connections between origin i and the first stop of each of the R(j) 
routes, all the connections between the pairs of bus stops, and all 
the connections between the last stop of each of the R(j) routes 
and the terminus are made in the least time consuming manner. A 
computational procedure for determining the lower bound upon the 
total traveling time of a set of routes is described in section V-B. 
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The c(i,j) are based upon a lower bound instead of an upper bound 
because an attempt is being made to develop a set of routes whose 
total traveling time approaches this lower bound. 
c(i,j) = 0 when school j is a fictitious school. 

The origin associated with the maximum x(i, j) for school j is 
defined to be the super-origin for school j. In case of a tie, the 
super-origin is selected arbitrarily. 

The allocation of a number of routes for school j to origin 1 
and the selection of a super-origin for each school serviced during 
the same period is determined once. 

C. Determination of Trial Routes: Nearest-City Approach 

A trial route which starts at the super-origin, visits every 
stop once, and terminates at the school, the route that an infinite 
capacity bus would traverse, is determined either by the Nearest City 
Approach or by Algorithm A discussed in the next section. The Nearest 
City Approach can be described as follows: 

Let n = number of stops assigned to the school (including 

origin and terminus) 
stop 1 = origin of route 
stop n = terminus or school 

r = number of the trial route being generated 
where r = 1,2 (n-2). 

The subscript denoting the particular school has been omitted in 
order to simplify the nomenclature. 

This trial route generated is one in which an "infinite capacity" 
bus starts at the super-origin, proceeds to stop (r+1) and then 
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repeatedly selects as its next stop that point which is nearest to its 
present stop and which has not yet been serviced until it reaches the 
terminus or school. The current best trial route, the one requiring 
the least traveling time, is determined and saved each time that this 
process is executed. Since this procedure specifies the first bus 
stop visited, a different trial route is generated each time. If the 
acceptance criterion has not been satisfied after (n-2) trial routes 
have been generated, then Algorithm A is used to develop succeeding 
trial routes from the best trial route previously determined, 

D. Determination of Trial Routes: Algorithm A 

Algorithm A, a systematic procedure for decreasing the total 

time required to traverse the infinite capacity bus route, is an 

29 30 

extension of Algorithm 1 developed by the author ’ . After the 

Nearest City Approach has been exhausted, the trial routes are 
generated by applying Algorithm A to the best infinite capacity bus 
route available. Initially Algorithm A is applied to the best 
route determined by the Nearest City Approach. Thereafter, Algorithm 
A is applied to the most recent trial route it generated. 

This algorithm determines sets of three links which can be 
changed simultaneously without destroying the continuity of the tour, 
the non-coincidence of the origin and the terminus of the route, and 
the direction of the unchanged portions of the route. The latter 
constraint is necessary, because the procedure is applicable to non- 
symmetric as well as symmetric problems. If the time required to 
traverse the three new links is less than the time required to 
traverse the links which they replace, then the new route becomes the 
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next trial route and the best infinite capacity bus route available. 
If the proposed change results in no improvement, then another set of 
changes is determined and examined. Algorithm A is repeated until it 
is unable to improve the best trial route available, i.e. Algorithm A 
is exhausted. 

The sets of three link changes in the current tour through a 
network of n points or stops in which the origin and terminus do not 
coincide are generated as follows for all combinations of i and j 
where 1 s i i (n-1), 2 jS j (n-1), stop 1 is the super-origin 

and stop n is the terminus or school: 

1. New link 1 starts at point i and ends at point j where 
i^j and point j ^ the point which follows point i in the 
current tour. 

2. New link 2 starts at point k and ends at the point which 
follows point i in the current route where kAi« 

k is cycled as follows: 
point k^ = point j 

point kj = point which follows point j in the 

• current route 

point k^ = point which lies (L-l) consecutive 

positions after point j in the current 
route in a clockwise direction 
point k^ = point which precedes point i in the 
current route 

In order . to determine the complete range of values for 
index k associated with a particular i,j combination, it 
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is necessary to assume the existence of a fictitious link 
between the terminus and origin of the current route while 
counting the consecutive positions after point j in a clock- 
wise direction. 

3. New link 3 starts at the point which precedes point j in the 
current route and ends at the point which follows point k in 
the current route. 

For illustrative purposes. Algorithm A will be applied to a route 
containing six points for one combinatior of i and j and the entire 
range of index k associated with it. Point one is the super-origin 
and point six is the terminus or school. The number associated with 
each node is the permanent identification number of the stop and 
corresponds to its position in the interstop travel time matrix, eg. 
stop four data would form row four of the interstop travel time 
matrix. 

Let i = 2 and j = 3 be the i,j combination 
Kq = 1-4- 2-5 -3- 6 be the current route 



4 2 
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The first value assumed by k is the number of point j or three. 
The proposed route is = 1-4-2-3-5-6. Links 2-5, 5-3, and 3-6 are 
replaced by new links 2-3, 3-5, and 5-6, indicated by the dashed arcs. 




The second value assumed by k is the number of the point which 
follows point j or six. However, no new route can be proposed by 
Algorithm A for this combination of i,j, and k because the value 
assigned to index k is the number associated with the terminus of the 
route. Since k has not yet assumed the value of the identification 
number of the point preceding point i in the current route, four, the 
cycle for index k is incomplete. 

The third value assumed by k is the number of the point which 
lies two consecutive positions after point j, in a clockwise 
direction, or one. The proposed route is = 1-5-4-2-3-6. Links 
1-4, 2-5, and 5-3 are replaced by new links 1-5, 5-4, and 2-3. 
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The fourth value of k is the number of the point which lies three 
consecutive positions, in a clockwise direction, after point j in the 
current route or four. This is the last k which can be generated for 
the combination i = 2 and j = 3 because the point numbered four 
precedes point i in the current route. The proposed route is 

= 1-4-5-2-3-6. Links 4-2, 2-5, and 5-3 are replaced by new links 
4-5, 5-2, and 2-3. 
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When Algorithm A cannot improve the best trial route available, 
it is said to be exhausted. This best trial route which Algorithm A 
cannot improve is considered to be the quasi-optimal modified 
traveling-salesman route and the last possible trial route. A 
computational procedure for Algorithm A is described in section 
V-A. 

E. Partitioning Procedure 

The partitioning procedure is applied to every trial infinite 
capacity bus route determined by either the Nearest City Approach or 
Algorithm A. This procedure requires the following additional 
information: 

1. a student load vector specifying the number of students 
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assigned to each stop, 

. the bus capacity, and 

3. the maximum allowable riding time of the students picked up 
at the first stop of any route. 

The partitioning procedure generates a set of bus routes each of 
which starts at the super-origin, visits the stops of the trial 
infinite capacity bus route in the order previously determined until 
the bus is loaded properly and proceeds to the terminus or school. 

At each stop, the bus load count is incremented by the appropriate 
element of the student load vector and the time tally is incremented 
by the traveling time from the previous stop. When either the bus 
capacity or the riding time constraint is about to be exceeded, the 
previous stop becomes the last one serviced by the bus before 
proceeding to the school. The next route starts at the super-origin 
and proceeds directly to that stop of the trial infinite capacity 
bus route which immediately follows the last stop serviced. All 
individual bus routes are determined in the same manner. The 
sequence of stops generated by either the Nearest City Approach or 
by Algorithm A is preserved throughout this procedure. If the set 
of routes determined by the partitioning procedure contains no more 
than the maximum allowable number of routes, then the improvement 
process, described in the next section, is applied to each route of 
the set. Otherwise, another trial infinite capacity bus route is 
generated by either the Nearest City Approach or Algorithm A. 

F. Improvement Process 

Each individual route belonging to a feasible set of routes 



developed by the partitioning procedure is then improved by application 
of Algorithm A, previously described in section IV-D. 

First, form a submatrix of the given interstop traveling time 
matrix consisting of those elements associated with the super-origin, 
the stops serviced by the individual route, and the terminus or school. 
The infinite capacity bus route corresponding to this submatrix is 
the individual route being improved. Then Algorithm A is applied to 
this route until it can make no further improvement. This two step 
procedure is repeated until all the individual bus routes of the 
feasible set have been improved. The given matrix of interstop travel 
times between each possible origin and every stop assigned to the 
school and the given matrix of interstop travel times for all pairs 
of stops except those links involving an origin are preserved at all 
times . 

G. Acceptance Criteria 

The best available set of routes, all of which start at the 
super-origin, for a particular school j will be considered to form 
the basis of the quasi-optimal solution if one of the following 
criteria is satisfied: 

1. The total time required to traverse this set of routes is 
less than or equal to the product of a factor specified by 
the transportation director, greater than one, and the 
lower bound on the total time required to traverse the 
number of routes contained in the best available set of 
routes, assuming that all of them start at the super-origin. 

. The total number of trial routes generated is about to 
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exceed some number specified by the transportation director. 
Thus, the school administrator has an opportunity to specify the 
degree of optimality required for the set of routes accepted as a 
final solution. In the event that neither of the criteria are 
satisfied and no further trial infinite capacity bus routes can be 
developed because both the Nearest City Approach and Algorithm A have 
been exhausted, then the best available set of routes will become the 
basis of the quasi-optimal solution by default. The latter situation 
is one in which the capability of the method is unable to satisfy the 
requirements of the school. 

H. Final Allocation of Origins to Individual Routes 

The feasible set of individual bus routes which satisfy the 
acceptance criterion specified by school j all start at the super- 
origin. If all the routes for school j should start at the super- 
origin, as previously determined, then a final allocation of origins 
to individual routes is unnecessary. However, if all the routes for 
school j do not start at the super-origin, then the remaining origins 
are allocated in such a way that the additional number of time units 
traveled will be minimum. This final allocation problem can also be 
expressed as a transportation problem. 

Let nl = number of origins supplying buses for the period 

n3(j) = number of routes belonging to the feasible set of 
routes satisfying the acceptance criterion 
specified by school j 

x(i, j) = number of routes which should start at origin i 



and end at school j 



d(i,j,k) = total traveling time of route k assigned to 
school j when it starts at origin i 
y(i,j,k) = ( 1 when route k assigned to school j starts 
"S Ht origin i 
\£> otherwise 

The problem can be stated as follows: 

Find variables y(i,j,k) for a fixed j and all combinations of i and k, 
where i = l,2,...,nl and k = 1,2, . . . ,n3(j) , such that the objective 
function 
nl n3(j) 

Y C d (i,j,k)] [y(i, j,k)] 
i=l k=l 

is minimized and restrictions one through four, stated below, are 
satisfied. Optimization is with respect to minimizing the total 
traveling time required by the n3(j) routes accepted by school j. 

The first constraint, 
nl 

Y y( i >J> k ) ■ for k = 1,2, . . . ,n3(j) , 

i=l 

insures that a route is assigned to only one origin. There will be 
n3(j) constraints of this type. 

The second constraint, 
n3(j) 

Y yUtJ.k) = x (i> j) > for i = 1,2, ...,nl 
k=l 

insures that the number of bus routes starting at origin i equals 
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the number of louses previously assigned to school j from origin i. 
There will be nl constraints of this type. 

The third constraint, 
y(i»j»*0 = 0 or 1, for all i and k, 
is self-explanatory. 

The fourth constraint, 
nl 

Y x (i> j) * n3(i) 

i=l 



insures that a feasible solution exists. Since the first constraint 
n3(j) nl n3(j) 

implies Y ^ y(i,j,k) = / 1 = n3(j) and the second 

k=l i=l k=l 

nl n3(j) nl 

constraint implies ^ ^ y(i,j,k) = Y x(i,j), then indeed 



i=l k^l 



i=l 



nl 



the fourth constraint must be satisfied. Since Y x(i, j) equals the 



i=l 



nl 



maximum allowable number of routes and n3(j) £ ^ x(i,j). 



i=L 



fictitious route may have to be introduced to use the extra available 
bus in order to satisfy the fourth constraint. 

The parameter nl was given; the parameter n3(j) was determined 
by the partitioning procedure described in section IV-E. The values 
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of x(i,j), the number of bus routes which should start at origin i 
and end at school j, were determined by solving the transportation 
problem described in section IV-B. The d(i,j,k) are calculated as 
follows: 

Let t(s,j,k) = total time required to traverse route k when it 

starts at the super-origin and ends at school j 
l(s,j,k) = traveling time between the super-origin and the 
first stop serviced by route k associated with 
school j 

l(i,j,k) = traveling time between origin i and the first 
stop serviced by route k associated with 
school j 

d(i,j,k) = t(s,j,k) - l(s,j,k) + l(i,j,k) 
d(i,j,k) = 0 when k is a fictitious route. 

After the final allocation of origins to the individual bus 
routes for school j has been completed, a lower bound upon the total 
traveling time of the number of routes required by school j is 
calculated. This procedure is described in section V-B. 
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V. COMPUTATIONAL PROCEDURE 



A. Algorithm A 

A computational scheme for Algorithm A is described by the flow 
chart given in Figure 1. The nomenclature used in the flow chart is: 

N « total number of bus stops (including the origin and 
the terminus) 

[m] = a square matrix of order (N+l) composed of an NxN 
interstop travel time matrix augmented by an 
additional row and an additional column 
M(I, j) *= number of time units required to travel from bus 

stop I to bus stop J: l^I^N.l^J^N, I^J 

M(I,I) = number of time units required to travel from bus 

stop I to the stop which immediately follows it in 
the current modified traveling-salesman route 
M(I,N+1) = identification number of the bus stop which im- 
mediately follows bus stop I in the current 
modified traveling- salesman route 
M(N+1, J) = identification number of the bus stop which im- 
mediately precedes bus stop J in the current 
modified traveling-salesman route 
M (N+l, N+l) is not used 

01 = origin of new link 1 
T1 = terminus of new link 1 
02~=~origin of new link 2 
T2 = terminus of new link 2 
03 = origin of new link 3 
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1 V 


01=1, TI=J 
02=J, T2=M(I,N+I) 
03=M(N+U), T3=M(J,Ntl) 
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T3 = terminus of new link 3 



The first row of [m] is related to the origin and row N is related to 
the terminus of the modified traveling-salesman route. 

Blocks 1-4 provide for the initialization and incrementation of 
index I and index J. 

Blocks 5-6 prevent the creation of an illegal combination of I 
and J. Block 5 prevents a change in which indexes I and J are identical. 
Block 6 prevents a change in which stop J immediately follows stop I in 
the current tour and thereby eliminates the possibility of generating 
a new route that is identical to the old one. 

Block 7 initializes the origin and terminus of each of the three 
new links associated with the current I,J combination. 

Block 8 calculates the time required to traverse the three old 
links which are candidates for replacement and the time required to 
traverse the proposed new links. 

Block 9 determines whether the proposed set of new links reduces 
total transit time. y 

Block 10 tests whether index k, the origin of new link 2, has ' 
assumed all possible values for the current combination of indexes I ^ 
and J. 

Block 11 tests whether the next value of index k would be the 
identification number assigned to the terminus of the current route. 

Block 12 determines whether it would be possible to assign another 

value to index k by introducing a fictitious link between the terminus 

... . , W 

and origin of the current route. 

Block 13 calculates the new value of index k, the origin of new 



er|c 



45 



link 2 as described in Section IV-D, and the terminus of new link 3 
when index k does not assume the value of the identification number of 
the terminus of the current route. 

Block 14 calculates the new value of index k, the origin of new 
link 2, and the terminus of new link 3 after a fictitious link was 
assumed to exist between the origin and terminus of the current route. 

Blocks 15-3.6 determine whether indexes I and J can be further 
updated. If I cannot be updated, then Algorithm A has been executed 
for all possible combinations of I and J. If Algorithm A is exhausted 
while it is being used to generate a trial infinite bus capacity route 
prior to execution of the partitioning procedure, then the best avail- 
able set of routes will have to be accepted as the quasi- optimal 
solution for the school. Exit will be to the final allocation of 
origins procedure. If Algorithm A is exhausted while it is being used 
to improve an individual bus route, then exit is to the improvement of 
the next route of the set. 

Block 17 incorporates the three new links, that have been found to 
reduce total transit time, into the current modified traveling-salesman 
route by altering row (N+l) and column (N+l) which store the sequence 
of stops in the new tour and by inserting the new transit times into 
the main diagonal elements. If Algorithm A produces an advantageous 
change while it is being used to create a new trial infinite capacity 
bus route, then exit is to the partitioning procedure. If Algorithm A 
is being used to improve an individual bus route when an advantageous 
change occurs, then Algorithm A is restarted and exit is to Block 1. 
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B. Lower Bound 



The lower bound on the total time required to traverse the set of 
routes accepted as the quasi-optimal solution by the school is calculated 
under the assumption that all the connections between an origin and the 
first stop of any route, all the connections between the pairs of bus 
stops, and all the connections between the last stop of any route and 
the terminus are made in the least time consuming manner. 



Let N = total number of bus stops (including origin and terminus) 

M = rectangular matrix with (N-l) rows and N columns 
M(I,J) = number of time units required to travel from bus stop 
I to bus stop J: 2^I^N, l^J^N, I^J 

Nl = number of origins supplying buses for the school 
Ml = rectangular matrix with Nl rows and N columns 
M1(K,J) = number of time units required to travel from origin K 
to bus stop J: 1 £ K £ Nl, l£j£N 
NR(K) = number of routes starting at origin K: 1 ^ K £ Nl 

NL 

Q = Y, MW 
K=1 

[M] is the given matrix of inteistop travel times for all pairs 
of stops except those links involving an origin. Since the origin of 
a route is not allowed tc be the terminus of a link, M(I,1) = “ for 

2 £ I s; N. Moreover, M(N,J) = <*> for 1 £ J £ N because the terminus of 

a route is not allowed to be the origin of a link. These two restric- 

tions are necessary because the method of solution is applicable to the 
general routing problem in which the origin and terminus do not coin- 
cide. 
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[Mil is the given matrix of interstop travel times between each 
origin servicing the period and every stop assigned to the school. 
Since no link is allowed between the origin and terminus of a route, 

Ml (K,N) ■ •tel £ K £ Nl. Moreover, M(I,I) = ® f or 2 £ I £ N and 
Ml(K,l) ■ ®for 1 £K £ Nl because no loop is permitted at any bus 
stop. 

Parameters Nl and NR(K) where 1 £ K £ Nl and NR(K) ^ 0 are deter- 
mined by the final allocation of origins to routes procedure described 
in Section IV-H. 

The lower bound is calculated as follows: 

1. Form [M2], a square matrix of order (Q+N-l), by the 
procedure described below. 

2. Reduce [M2 ] until there is at least one zero in every 
row and column. This is accomplished by subtracting the 
smallest element in each row from every element in the 
row, and then subtracting the smallest element in each 
column of the remaining matrix from every element in the 
column. The lower bound on the final set of routes is 
the total reduction or the sum of the elements subtracted 
from the rows and columns. 

Starting at row one, [M2] is formed as follows: 

1. For each K, 1S R S Nl, row K of [Ml] is stored NR (K) times 
in [m 2], i.e. each route starting at origin K contributes 
one row to [m 2],. This rectangular submatrix consisting of 
Q rows and N columns occupies the upper left portion of 
[M2]. During the reduction process, this submatrix provides 
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that portion of the lower bound contributed by traveling 
between any origin and the first stop of a route. 

2. [m] is stored in rows (Q+l) through (Q+N-l) of [M2], 

During the reduction process the first (N-l) columns 

of this submatrix provide that part of the lower bound 
contributed by traveling between any pair of points. 

3. Steps one and two create columns one through N of [m2]. 

Column N is then stored (Q-l) more times in [m 2], i,e. 
columns N through (Q+N-l) of [M2 ] are the same vector 
whose elements represent the time required to travel 
between any bus stop and the terminus. During the re- 
duction process, the aubmatrix occupying rows one through 
(Q+N-l) and columns N through (Q+N-l) of [m 2] provides 
that portion of the lower bound contributed by traveling 
between the last stop of any route and the terminus. 

When the lower bound has to be calculated for a set of routes, all 
of which start at the same origin, then Ni = 1 and NR(1) = Q equals the 
number of routes contained in tne set for this procedure. 
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VI. EVALUATION 



Any heuristic procedure should be evaluated on two bases. First, 
in order to be acceptable, the method must be reasonable with respect 
to the size of the problem it can handle, the assumptions it imposes 
and the logical processes it uses. Second, in order to be practical, 
the procedure must be able to yield answers to the problem being 
considered at a cost commensurate with the value received. 

A. Model and Method 

The model developed in this dissertation is a general one which 
utilizes variables that are applicable to all school systems. Moreover, 
it requires few assumptions- that simplify the real world environment. 

The method used to obtain a solution for the model reduces the routing 
of buses for an entire school district to a set of sequential steps that 
can be readily programmed for a digital computer. Furthermore, this 
set of sequential steps is arranged into groups called iterations or 
passes such that a "reasonably good" feasible routing system can be ob- 
tained site r a few iterations have been completed. Thus, the user of 
the model and method can specify the degree of optimality desired for 
the routing system being developed. 

The procedure used to develop a set of bus routes for each school 
is logical. At the beginning of a period, the available buses at the 
origins are allocated to every school serviced during the period so that 
the total estimated traveling time required by all the routes developed 
for the period is minimized. This is accomplished by solving a trans- 
portation problem in which the elements of the requirements vector 
are the maximum allowable number of buses required by a school, the 
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elements of the availability vector are the number of buses located 
at each origin at the beginning of the period, and the elements of the 
cost matrix are based upon the lower bound on the total traveling time 
required by each combination of origin and school. Since the bus 
capacity and riding time constraints often make the generation of a 
feasible set of routes containing the absolute minimum number of routes 
unattainable, the selection of this requirement vector is considered 
reasonable. Moreover, because an attempt is being made to develop a 
set of routes whose total traveling time approaches the lower bound, the 
estimation of the elements of the cost matrix used by the method is a 
logical choice. 

The determination of a trial route which starts at the "super- 
origin", visits every stop once and terminates at the school by either 
the Nearest City Approach or Algorithm A tends to group stops that are 
located in the same neighborhood. Since a school usually services an 
area of less than twenty square miles, the variability associated with 
the magnitude of the elements of the inter3top travel time matrix is 
low. Thus, arranging the bus stops to be serviced into groups located 
in the same neighborhood is logical. The first (n-2) different trial 
routes, where n is the total number of stops assigned to the school 
including the origin and the terminus, are developed by the Newest 
City Approach because it is a rapid and relatively efficient process. 
After all trial routes are generated by the Nearest City Approach, the 
route requiring the least traveling time is saved in order to reduce 
the number o£ trial routes which Algorithm A will generate before it 
is exhausted and to insure that Algorithm A does not create one of the 
routes previously developed. 
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Algorithm A is applied to the best available trial route because 
a good "infinite capacity'.' bus route usually partitions into a good 
routing system even though the best "infinite capacity" bus route does 
not necessarily generate the best set of individual routes. Although 
Algorithm A is a slightly slower process than the Nearest City Approach, 
it guarantees that the next trial route developed will require less • 
traveling time. Both of the algorithms were selected because they are 
rapid, efficient, and exhausted in a finite number of steps. 

The partitioning procedure quickly generates routes which not only 
satisfy the bus capacity and passenger riding time constraints but also 
service stops in the same general area, if the set of routes developed 
contains at most one more than the absolute minimum number of routes 
required by the school, then each route of the set is improved by ap- 
plication of Algorithm A until it is exhausted while preserving the 
assignment of stops to the route made by the partitioning process. This 
improvement procedure may reduce the traveling time required by the 
individual route and thus tends to reduce the total transit time of the 
routing system. 

Minimizing the number of routes included in a routing system in- 

\ 

sures that the buses are being utilized to full capacity and also tends 
to reduce the total traveling time required by the set of routes by 
eliminating some of the links between an origin and the first stop of 
any route and the links between the last stop of any route and the 
terminus. Since at most ten students are usually assigned to a bus 
stop in order to reduce the noise level and the possibility of landscape 
damage in residential areas, satisfying the maximum allowable number of 
routes in a routing system constraint presents no problem. However, if 
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a transportation director decides to lump the bus stops extensively, so 
as to reduce the size • of the interstop travel time matrix or if the 
population density is such that many students are assigned to one stop, 
then determining a feasible routing system containing at most one more 
than the absolute minimum number of routes may be unattainable. The 
user of this method must then adjust the set of constraints by either 
increasing the maximum allowable number of routes in the routing system 
or by reducing the bus capacity. Another alternative is the addition 
of bus stops to reduce the number, of students assigned to individual 
stops. 



The specification of a criterion for accepting a set of routes 
gives the transportation director the opportunity to select the degree 
of optimality under which the routing system will be developed. This 
acceptance criterion is applied to a routing system when all the routes 
start at the "super -origin" because allocating an origin to each route 
every time a feasible routing system is generated would require more 
computer time than is warranted by the improvement which would be 
realized. 

The final allocation of origins to individual routes, if necessary, 
is accomplished by solving a transportation problem in which the elements 
of the availability vector are the number of routes which start at 
each origin as determined by step one of the solution procedure and the 
elements of die requirements vector are all equal to one. Each element 
of the cost matrix is the total traveling time required by a route when 
it starts at a particular origin. Optimization is with respect to 
minimizing the total traveling time required by the final set of routes 
for the school. Initially, the solution procedure develops routing 
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systems in which all the individual routes start at the "super-origin" 
in order to minimize the effects of this final allocation of origins to 
routes. After an origin has been assigned to each route of the set, no 
attempt is made to further reduce the transit time required by each 
route because the improvement process may create a violation of the 
student riding time constraint. 

Although the final lower bound on the set of routes considered to 
be the quasi-optimal routing system does not include the effects of the 
bus capacity and riding time constraints, it does include a partial 

J ‘ 

effect of the number of routes included in the system. However, it 
provides some measure for assessing the quality of the set of routes 
developed even though it may be unattainable. 

Thus, this heuristic procedure is considered to have a sound basis* 

B. Computational Experience 

Since there is no known bus routing method suitable for use on a 
computer which will guarantee an optimum solution, any heuristic pro- 
cedure must be judged not only with respect to its degree of success 
relative to the best set of routes available for known problems but also 
with respect to its consistency in the level of success. 

A computer program based upon the procedure described in this 
dissertation was written in FORTRAN IV for use on the CDC 6400 computer-. 
It can handle bus routing for a school district involving any number of 
periods, nine route origins and nine schools per period*’ one hundred 
twenty stops per school, thirty-four routes per school, and thirty 
stops per route. All computations are done in integer arithmetics 

Computational experience was gained in three phases. First, 
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Newton's procedure was applied to school bus routing and delivery pro- 
blems reported in the literature. Then, the relationship between the 
amount of computer time required to develop a routing system and the 
number of stops serviced by the set of routes was examined. Finally, 
this method was used to develop a set of bus routes for four schools 
in the Williamsville Central School District, a suburban area in 

Western New York. 

34 

Thompson drew a map of a hypothetical school system showing the 
location of the school, the location of 31 bus stops specified by the 
school, and the number of students assigned to each stop. Moreover, 
all roads were marked off in units of one-half mile. In order to 
simulate the area of a real school system as closely as possible, the 
map included features such as: isolated areas, contour roads, and 

varying population densities.' The elements of the symmetric interstop 
distance matrix calculated from this map ranged between 1 and 13.5 miles; 
the elements of the student load vector varied between 3 and 22 students. 
Although it was stated that 3 minutes were required to travel one mile 
and that loading at each bus stop required a minute, no restriction was 
imposed on the student riding time. Buses of 30, 36, 42, 48, 54, 60, 

66, and 72 passenger capacities were available. The problem was to 
design the set of bus routes, all of which started and ended at the 
school, required to transport the 252 students assigned to the 31 bus 
stops. 

The group of 49 school superintendents with experience in school 
bus routing and the group of 37 transportation directors who participated 
in Thompson's study were allowed to design the routing system for the 
hypothetical school by any manual method they favored. 



The best set of routes developed by a member of the group of tran- 
sportation directors involved 5 routes whose loads varied between 28 
and 67 students and whose total traveling distance was 100 miles . The 
best set of routes designed by a member of the group of school super- 
intendents consisted of 5 routes whose loads ranged between 33 and 67 
students and whose total traveling distance was 101 miles. 

Newton's method was applied to the same hypothetical school for all 
combinations of 4 bus capacities and 4 maximum allowable riding distances 
or a total of 16 cases. Each case was run until Algorithm A was ex- 
hausted at iteration or pass 49, i.e. the Nearest City Approach generated 
31 trial routes and Algorithm A generated 17 "infinite capacity" bus 
routes. 177.133 seconds of computer time were used to develop the rout- 
ing systems for the 16 c^ses or approximately 11 seconds per case. Most 
of the solutions were accepted from one of the last 3 trial routes ge- 
nerated by Algorithm A as was expected. A summary of the results 
appears in Table 1. 

Case 13 defined by a bus capacity of 54 students and a maximum 

riding distance of 20 miles yielded the best routing system whose total 

traveling distance was 94.5 miles. This represents a savings of 5.8% 

miles with respect to the best set of routes developed by a transportation 

34 

director participating in Thompson's study . The best set of routes 
for Case 13 is: 
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THOMPSON DATA 



31 STOPS SERVICED 

252 STUDENTS TRANSPORTED 

APPROXIMATELY 11 SECONDS CDC 6400 COMPUTER TIME PER CASE 



CASE 


CAPACITY 

OF 

BUS 


MAXIMUM 

RIDING 

DISTANCE 


MINIMUM 

NUMBER 

ROUTES 


NUMBER 

ROUTES 

USED 


DISTANCE 

ROUTE 

SYSTEM 


RATIO 


NUMBER 

ANSWER 

PASS 


1 


72 


20.0 ml. 


4 


5 


97.0 ml. 


1.48 


44 


2 


72 


18.0 


4 


5 


97.0 


1.48 


47 


3 


72 


17.5 


4 


5 


95.5 


1.46 


47 


4 : 


72 


17.0 


4 


5 


95.5 


1.46 


48 


5 


66 


20.0 


4 


5 


97.0 


1.48 


44 


6 


66 


18.0 


4 


5 


97.0 


1.48 


47 


7 


66 


17.5 


4 


5 


95.5 


1.46 


47 


8 


66 


17.0 


4 


5 


95.5 


1.46 


48 


9 


60 


20.0 


5 


5 


95.5 


1.46 


44 


10 


60 


18.0 


5 


5 


95.5 


1.46 


47 


11 


60 


17.5 


5 


6 


103.5 


1.54 


47 


12 


60 


17.0 


5 


6 


106.0 


1.58 


46 


13 


54 


20.0 


5 


5 


94.5 


1.44 


44 


14 


54 


18.0 


5 


6 


102.0 


1.52 


47 


15 


54 


17.5 


5 


6 


102,0 


1.52 


47 


16 


54 


17.0 


5 


6 


105.5 


1.57 
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Table 1 




l 




Number 



Route 



Load 



Miles 



1 


School-O-U-P-R-Q-RR-S-School 


50 


17.5 


2 


School- V-T-Z-W- School 


50 


16.0 


3 


School -Y-F-E-A-B-D-C-DD- School 


50 


23.0 


4 


School-X-G-H-HH- 1- J- School 


53 


21.5 


5 


School-OO-K-LL-M-L-N-School 


49 


16.5 



Cases 3, 4, 7, 8, 9 and 10 produced routing systems whose total 
traveling distance was 95.5 miles. However, all the routing systems were 
not identical. Cases 3, 4 , and 7 resulted in the same routing system; 
case 8 yielded a second set of routes; cases 9 and 10 produced a third 
routing system. Thus, by developing routing systems for various com- 
binations of bus capacity and passenger riding distance constraints and 
a fixed interstop distance matrix it may be possible to produce alter- 
native routing systems whose total traveling distance is identical. A 
routing system can then be selected from these sets upon the basis of 
either bus load or individual route length variability or some other 
statistic considered important by the school transportation personnel. 
Studying alternate optimal solutions can be easily accomplished by us- 
ing this computational procedure because of its speed and efficiency. 

Boyer ^ designed a set of bus routes, all of which start and end at 
the school, to transport 575 students assigned to 45 bus stops for a 
school in Hennepin County, Minnesota. The problem involved a symmetric 
interstop travel time matrix, a bus capacity constraint of 65 students, 
and a riding time constraint of infinity. Since Boyer's method does 
not require knowledge of all the elements of the interstop travel time 
matrix, only 171 elements were listed in (6). The elements of the 
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student load vector ranged between 1 and 33 students; the known elements 

f 

of the inters top travel time matrix varied between 1 and 23 minutes. 

. \ 

Although the absolute minimum number of buses needed to provide tran- 
sportation for the school was 9, Boyer developed a routing system using 
11 buses and! requiring 413 minutes traveling time. 

Newton's procedure was applied to the same problems until Algorithm 
A was exhausted. At iteration 51, i.e. the Nearest City Approach ge- 
nerated 45 trial routes and Algorithm A generated 5 trial routes. A 
routing system using 10 buses whose total traveling time was 394 minutes 
with a lower bound of 221 minutes was developed in 13.392 seconds. The 
quasi- optimal routing system, obtained from iteration 49 of this pro- 
cedure, requires 4.8% fewer minutes than Boyer's solution. The author 
feels that a better solution would have been obtained with the complete 
interstop travel time matrix. The best set of routes is: 



Number 


Route 


Load 


Time 


1 


S-22-23-S 


60 


32 


2 


S-24-25-26-31-21-20-19-S 


65 


38 


3 


S-37-41-36-35-S 


65 


27 


4 


S-40-39-34-38-S 


60 


30 


5 


S-33-32-29-30-27-28-12-S 


57 


62 


6 


S-6- 13-11-5- 1-S 


63 


49 


7 


S-2-3-4-7-9-8-S 


56 


48 


8 


S- 10-43-42-S 


62 


36 


9 


S-44-45-18-17-S 


61 


36 


10 


S-16-15-14-S 


26 


36 



Thus, Newton's method was able to obtain better routing systems 
for the only two school bus routing problems which have appeared in the 
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34 

literature. The results for Thompson's data -3 were especially encourag- 
ing because experienced transportation directors can usually develop 
almost optimum routing systems for problems involving 20-30 stops by 
visual trial-and-error adjustment. 

The rest of the problems from the literature which were solved by 
this computational procedure are delivery problems. Although, the 
school bus scheduling problem and the delivery problem are conceptually 
the same, differences which exist in them must be considered when 
evaluating a heuristic procedure. 

The school bus scheduling problem is characterized by a large 
number of stops which must be serviced and a non-symmetric interstop 
travel time matrix for reasons such as: restrictions imposed upon the 

crossing of busy streets by students, one-way streets and limited access 
highways. The elements of the interstop time matrix usually have a 
narrow range and a low variability because of the relatively small area 
serviced by a school and the restrictions imposed upon student walking 
time to the bus stop. The elements of the student load vector also have 
a narrow range and a low variability because of upper bounds usually 
placed upon the number of students assigned to a bus stop to reduce the 
noise level and the possibility of landscape damage in residential 
areas. Moreover, the bus stops are usually arranged in groups because 
of housing developments and the placement of bus stops along main 
thoroughfares in sparsely populated areas. Even in sparsely populated 
areas, no bus stop is really isolated, i.e. each stop is a relatively 
short distance from either the school or any other stop designated by 
the school. Moreover, any computer based procedure for the school bus 
scheduling problem must require little computer time in order to be of 



practical value to a school district because of the frequent updating 
of routing systems necessitated by population fluctuations, changes in 
school boundaries, and the building of new schools and lack of funds. 

The delivery problem is characterized by relatively few stops and 
a symmetric interstop distance matrix whose elements usually have a 
wide range and high variability because of the large area serviced by a 
warehouse or depot. Since no restrictions are placed upon customer 
demands, the elements of the customer demand vector also have a wide 
range and high variability. Moreover, the delivery stops are not 
necessarily arranged in groups because customer demand is not area de- 
pendent. Then too, any computer based procedure for the delivery 
problem may require a great deal of computer time and still be accept- 
able to a corporation which normally allocates ample funds for the 
development of routing systems. 

These differences between the delivery problem and the school bus 
scheduling problem are great enough to make the widespread interchange 
of heuristic solution procedures infeasible. Therefore, an efficient 
computational procedure which was designed primarily to handle routing 
problems possessing the characteristics of the school bus scheduling 
problem is not expected to be consistently superior with respect to the 
quality of the routing systems developed when applied to delivery 
problems. 

14 

Dantzig and Ramser developed a set of truck routes, all of which 
start and end at the depot, to deliver 18200 gallons of material to 12 
customers in 6000 gallon capacity trucks. The elements of the symmetric 
interstop distance matrix varied between 5 and 52 units; the elements 
of the customer demand vector ranged between 1100 and 1900 gallons. 
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The Dantzig and Ramser method produced a routing system using 4 routes 

whose total traveling distance was 294 units. 

9 

Clarke and Wright developed a routing system for the same problem 
requiring 4 routes covering 290 units, the conjectured optimum. 

Newton's procedure was applied to the same delivery problem until 
Algorithm A was exhausted at interation 13. A routing system using 4 
routes whose total traveling distance was 304 distance units was de- 
veloped in 1.117 seconds. The solution, obtained from iteration 5 of 
this procedure, requires 3.4% more distance units than the Dantzig and 
Ramser solution and 4.8% more distance units than the conjectured 
optimum. 

The set of routes developed by this procedure is: 



Number 


Route 


Load 


Distance 


1 


0-6-7-5-0 


4300 


64 


2 


0-9-8-10-0 


5300 


92 


3 


0-11-12-4-3-0 


5700 


120 


4 


0-2-1-0 


2900 


28 



9 

Clarke and Wright designed a set of truck routes, all of which 
start and end at the depot, to deliver 104,300 pounds of goods to 30 
customers in 14,000 pound capacity trucks. The elements of the 
symmetric interstop distance matrix ranged between 3 and 98 miles; the 
elements of the customer demand vector varied between 100 and 12,300 

9 

pounds. The Clarke and Wright solution used 8 trucks, the absolute 

minimum number of trucks possible, and required 1427 miles. For the 

lA 

same problem, the Dantzig and Ramser method produced a routing system 
using 10 routes whose total traveling distance was 1766 miles. By 
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18 

visual trial-and-error adjustment, Gaskell , developed a routing system 
for this problem using 8 routes and covering 1416 miles. 

Newton’s procedure was applied to the same problem until Algorithm 
A was exhausted at iteration 44. A routing system using 9 routes whose 
total traveling time was 1544 miles was developed in 7.948 seconds. 

The solution, obtained from iteration 5, required 9% more miles than 
the Gaskell solution, 8.2% more miles than the Clarke and Wright method, 
and 14.3% fewer miles than the Dantzig and Ramser solution. 

The set of routes developed by this method is: 



Number 


Route 


Load 


Distance 


1 


0-6-5-11-16-15-9-7-13-29-0 


13700 


215 


2 


0- 12- 14-4-3-24*.25U.23-0 


13500 


207 


3 


0-27=26-0 


1 1 onn 


199 


4 


0-8-10-19-0 


8700 


176 


5 


0-18-25-20-0 


12200 


150 


6 


C- 2-1-21-17-0 


8500 


111 


7 


0-30-0 


12300 


136 


8 


0-28-0 


9500 


186 


9 


0-19-0 


14000 


164 



Cochran^ designed the routing systems for two delivery problems 
using a modified Clarke and Wright method. Problem 1 involved 14,461 
units to be delivered to 12. customers in trucks of 4500 unit capacity. 
The elements of the symmetric interstop distance matrix varied between 
8 and 315 units; the elements of the customer demand vector ranged bet- 
ween 100 and 3726 units. The Cochran solution^ used 4 routes whose 
total traveling distance was 1433 units. 

Newton's method was applied to Problem 1 until Algorithm A was 
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exhausted at iteration 19. A routing system containing 4 routes whose 
total traveling distance was 1383 units was developed in 1.717 seconds. 
The solution, obtained from iteration 18 of this procedure, required 
3.6% fewer distance units than the solution produced by the modified 
Clarke and Wright method. Although the relative locations of the 
delivery stops are unavailable, the author feels that this computational 
procedure produced a superior solution because the delivery stops were 
arranged in groups and thus Problem 1 resembled a school bus scheduling 
problem. 

The set of routes developed by this method is: 



Number 


Route 


Load 

■' •.*> illinium* 


Distance 


1 


0-7-8-5-0 


3290 


185 


2 


0-14-0 


3726 


444 


3 


0-6-12-13-11-10-0 


3745 


478 


4 


0-4-3-2-9-0 


3700 


276 



The first two routes appeared in both routing systems. 

Cochran's Problem 2^ involved 1405 units of goods to be delivered 
to 25 customers in 120 unit capacity trucks. The elements of the 
symmetric interstop distance matrix varied between 2 and 221 distance 
units; the elements of the customer demand vector ranged between 15 
and 100 units. Using a modified Clarke and Wright method, Cochran 
designed a routing system using 14 trucks whose total traveling distance 



was 1468 distance units. 

This computational procedure was applied to Problem 2 until 
Algorithm A was exhausted at iteration 31. A routing system using 



trucks whose total traveling distance was 1486 distance units was 
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developed in 3.637 seconds. The solution, obtained from iteration 12 
of this procedure, travels 1.2% more distance units than the solution 
obtained by the modified Clarke and Wtight method. Both solutions used 
2 more trucks than the absolute minimum number of trucks possible. This 
was due to the combination of customer demand loads, most of which were 
greater than 50 demand units, and the location of the customers, i.e. 
customers in the same general area had total demands which exceeded the 
truck capacity and thus one route could service only one or two stops. 

The set of routes developed by this method is: 



Number 


Route 


Load 


Distance 


1 


0-13-11-0 


90 


82 


2 


0-7-15-0 


120 


70 


3 


0-16-14-0 


120 


102 


4 


0-2-0 


60 


4 


5 


0-3-0 


80 


10 


6 


0-4- 9-0 


110 


56 


7 


0-5-0 


90 


28 


8 


0-8-6-0 


115 


76 


9 


0-10-0 


60 


48 


10 


0-12-0 


90 


54 


11 


0-17-18-0 


120 


185 


12 


0-19-20-21-0 


110 


210 


13 


0-24-25-0 


120 


276 


14 


0-23-22-26-0 


120 


285 



Routes 1, 5, 6, 7, 10, 11, 12, 13, and 14 appeared in both routing 
systems. If routes 4 and 9 are merged into one route 0-2-10-0 which 
requires 52 distance units, then the routing system will use one less 
truck. However, the total traveling distance of the routing system with 
the merged routes in this case happens to remain unchanged. Although 



the purpose of a computer based routing procedure is to avoid visual 
trial-and-error adjustments which become less efficient and highly im- 
practical as problem size increases, this merging of two routes is 

noted as a point of interest. 

18 

Gaskell created four new delivery problems which were constrained 
by both truck capacity and route length. Total route distance included 
not only the number of units covered in traveling between delivery 
stops but also an allowance of 10 miles for each customer serviced. 

These problems were designed to compare the efficiency of variations of 
the Clarke and Wright method with respect to groupings of stops and 
isolated customers. 

Routing systems for each of these delivery problems were developed 
by the Visual Method, a combination of trial-and-error adjustment and 
manual permutation of groups of delivery stops. The routing systems 
produced by the Visual Method were the best available and therefore 
were used to evaluate the routing systems developed by the computer 
based variations of the Clarke and Wright method. Since these problems 
involved at most 32 delivery stops, the time consuming Visual Method 
produced optimal routing systems. However, for larger problems it would 
be unable to compete with computer based procedures. Thus, the routing 
systems developed by the computational procedure of this dissertation 
will be compared primarily with the routing systems produced by 
variations of the Clarke and Wright method, a computer based procedure. 

The first problem, Case Study Number 3, involved 29,370 units to 
be delivered to 32 customers in 8000 unit trucks. Bach route was not 
to exceed 240 miles including the mileage allowance of 10 miles per 
customer. The elements of the interstop distance matrix varied between 



1 and 118 miles; the elements of the customer demand vector ranged 
between 40 and 4000 units. This case is characterized by a close group- 
ing of some of the customers, a centrally located depot with respect 
to 30 of the delivery stops and 2 customers located a great distance 
from the depot. 

The Visual Method yielded a routing system using 4 routes whose 
total traveling distance was 813 miles. The best solution produced by 
a variation of the Clarke and Wright method used 5 rout* 3 and covered 
821 miles. However, the poorest solution produced by a variation of 
the Clarke and Wright method used 5 routes and required 850 miles. 

Newton ! s procedure was applied to Case Study Number 3 until 
Algorithm A was exhausted at iteration 39. It produced a routing system 
using 5 routes whose total traveling distance was 886 miles. The so- 
lution, obtained from iteration 23 of the procedure, required 7.9% more 
miles than the best solution obtained by a variation of the Clarke and 
Wright method. The reason that this computational procedure designed 
a less desirable routing system was the presence of the isolated de- 
livery stops which are generally not present in the school bus schedul- 
ing problem. 

The set of routes developed by this method is: 



Number 


Route 


Load 


Distance 


1 


0-17-24-23-22-20-21-18-19-15-14-0 


6900 


228 


2 


0-1-11-5-6-7-8-9-10-32-13-0 


7920 


177 


3 


0-31-30-3-4-2-12-0 


6350 


222 


4 


0-29-28-27-26-25-0 


7500 


147 


5 


0-16-0 


700 


112 




The second problem, Case Study Number 4, involved 22,500 units to 
be delivered to 21 customers in 6000 unit trucks. Each route was to be 
less than 200 miles in length including the mileage allowance of 10 
miles per customer. The elements of the interstop distance matrix 
varied between 3 and 83 miles; the elements of the customer demand vector 
ranged between 100 and 2500 units. None of the customers are isolated 
from either the depot or from other customers. The variability of the 
elements of the interstop distance matrix is relatively low. 

The Visual Method produced a routing system using 4 routes whose 
total traveling distance was 585 miles. The best solution produced by 
a variation of the Clarke and Wright method used 4 routes and covered 
598 miles. However, the poorest solution produced by a variation of the 
Clarke and Wright method used 4. routes and required 648 miles. 

Newton's procedure was applied to Case Study Number 4 until Al- 
gorithm A was exhausted at iteration 25. A routing system using 4 routes 
whose total traveling distance was 593 miles was developed. The solu- 
tion, obtained from iteration 22 of this procedure, required 1% fewer 
miles than the best and 9% fewer miles than the poorest routing systems 
designed by variations of the Clarke and Wright method. 

The set of routes developed by this method is: :i 



Number 



Route 



Load 



Distance 



1 

2 

3 0 : 

4 iK 



0 - 6 - 1 - 2 - 5 - 7 - 9-0 
0-10-8-3-4-11-13-0 
0-12-15-18-16-14-0 
0-17-20-21-19-0 



5600 

5400 

5500 

6000 



173 

162 

126 

132 



1 and 2 also appear in the routing system developed by the Visual 
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Routes 



Method. Since none of the customers were isolated and the variability 
of the elements of the interstop distance was low. Case Study Humber 4 
possessed two of the characteristics of the school bus scheduling pro- 
blem and this computational procedure developed a superior routing 
system. 

The third problem. Case Study Number 5, involved 12,750 units of 
goods to be delivered to 29 customers in 4500 unit capacity trucks. 

Each route was not to exceed 240 miles including the mileage allowance 
of 10 miles per delivery stop. The elements of the interstop distance 
matrix ranged between 1 and 121 miles; the elements of the customer 
demand vector varied between 100 and 3100 units. This problem was 
characterized by a loose grouping of customers located at various dis- 
tances from the depot. Gaskell considered this to be a difficult 
problem. 

The Visual Method designed a routing system containing 4 routes 
whose total traveling distance was 876 miles. The best solution pro- 
duced by a variation of the Clarke and Wright method used 5 routes and 
required 943 miles. However, the poorest routing system developed by a 
variation of the Clarke and Wright method involved 5 routes covering 
1017 miles. 

Newton's method was applied to Case Study Number 5 until Algorithm 
A was exhausted at iteration 33. A routing system using 5 routes whose 
total traveling distance was 913 miles was developed. The solution, 
obtained from iteration 30 of this procedure, required 3.3% fewer miles 
than the .best and. 11.3% fewer miles than the poorest routing systems 
designed by variations of the Clarke and Wright method. 
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The set of routes developed by this method is: 



Number 


Route 


Load 


Distance 


1 


0-21-14-8-9-17-12-11-10-28-18-0 


4125 


220 


2 


0-15-16-7-13-0 


1000 


156 


3 


0- 26-28-27-25-24-29-0 


2850 


234 


4 


0-3-6-1-4-5-2-0 


3975 


216 


5 


0-22-20-19-0 


800 


87 



Again, this computational method produced a superior routing system 
because the customers were arranged loosely in groups, a characteristic 
of the school bus scheduling problem. 

The fourth problem. Case Study Number 6, involved 10,189 units of 
goods to be delivered to 22 customers in 4500 unit capacity trucks. 

Each route was not to exceed a length of 240 miles including the mileage 
allowance of 10 miles per delivery. The elements of the interstop dis- 
tance matrix ranged between 4 and 145 miles; the elements of the 
customer demand vector varied between 60 and 4100 units. Some of the 
customers are loosely arranged in groups with the distance between 
customers greater than the distance between neighboring stops in other 
problems. One customer is isolated at a relatively great distance from 
the depot. 

The Visual Method produced a routing system involving 5 routes whose 
total traveling distance was 949 miles. The best solution produced by a 
variation of the Clarke and Wright method used 5 routes and required 
955 miles. However, the poorest solution produced by a variation of 
the Clarke and Wright method required 6 routes covering 1015 miles. 

Newton's procedure was applied to Case Study Number 6 until Al- 
gorithm A was exhausted at iteration 25. A routing system involving 6 



routes whose total traveling time was 1009 miles was developed. The 
solution, obtained from iteration 14 of this procedure, requires 5.7% 
more miles than the best solution produced by a variation of the Clarke 
and Wright method. 

The set of routes developed by this method is: 



Number 


Route 


Load 


Distance 


1 


0-14-17-15-16-3-2-0 


1144 


227 


2 


0-11-13-6-1-0 


775 


156 


3 


0-10-0 


4100 


78 


4 


0-12-9-5-4-8-7-0 


2700 


200 


5 


0-18-19-22-20-0 


1295 


216 


6 


0-21-0 


175 


132 



Routes 2, 3, and 6 appeared in both this routing system and the 
one developed by the Visual Method. 17.406 seconds of computer time 
were used to solve the last four problems. Again, as a point of in- 
terest, if routes 2 and 6 are merged to form the route 0-21-11-13-6-1-0 
which requires 275 miles, then the routing system will contain one less 
route and the total traveling distance will be reduced to 996 miles. 

Newton's procedure, developed primarily to handle the school bus 
scheduling problem, behaved as was to be expected when applied to eight 
delivery problems from the literature. When the delivery problem pos- 



sessed some of the characteristics of the school bus scheduling problem, 
this computational method produced a superior routing system. On the 
other hand, when the delivery problem involved isolated customers and 
many customers who were not even loosely grouped in areas, then this 
computational method produced a less desirable routing system than the 
other, computer based methods. However, the number of extra miles 
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required by the routing systems developed by this computational procedure 

never exceeded the best solution generated by a variation o£ the Clarke 

and Wright method by more than 8.2% which occurred in the Clarke and 

q 

Wright 30 stop problem . Thus, from a practical standpoint, this com- 
putational procedure performed at an acceptable level even when applied 
to delivery problems possessing characteristics which the method was not 
designed to handle. 

Since the ratio of the total traveling time or distance required by 
the set of routes developed divided by the lower bound on the total 
traveling time or distance is used not only to assess the efficiency of 
the routing system but to form a criterion for the acceptance or rejec- 
tion of a set of routes as the quasi-optimal solution for the problem 
being considered, an examination of this ratio with respect to the 
problems from the literature is warranted. 

A summary of the results obtained by applying this computational 
procedure to two school bus scheduling and eight delivery problems 
from the literature is given in Table 2. The column labeled % gives 
the percentage by which the method of this dissertation either exceeded 
or improved the total number of time/distance units required by the 
best routing system developed by a computer based procedure for each 

34 

problem being considered. In the case of the Thompson-31 stop problem , 
routing systems developed by manual methods were the only ones 
available. 

By observation, the ratio of the total traveling time/distance for 
a routing system divided by the lower bound on the total traveling 
time/distance yields little information about the relative efficiency 
of a routing system developed by this computational procedure in 
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SUMMARY OF RESULTS FOR PROBLEMS IN THE LITERATURE 



PROBLEM NAME 


TOTAL TIME/DISTANCE 
ROUTE SYSTEM 


ROUTE SYSTEM 
LOWER BOUND 


RATIO 


% 


j.THOMPSON-31 STOPS 


94.5 


65.5 


1.44 


-5.8 


BOYER-45 STOPS 


394,0 


221.0 


1.78 


-4.8 


JDANTZIG & RAMSER-12 STOPS 


304.0 


145.0 


2.10 


+4.8 


CLARKE & WRIGHT- 30 STOPS 


1544.0 


868.0 


1.78 


+8.2 


! COCHRAN- 12 STOPS 

.1 


1383.0 


816.0 


1.69 


-3.6 


COCHRAN-25 STOPS 


1486.0 


383.0 


3.88 


+1.2 


GASKELL #3-32 STOPS 

;'C .(;■ 


886.0 


620.0 


1.43 


+7.9 


GASKELL #4-21 STOPS 


593.0 


432.0 


1.38 


-1.0 


GASKELL #5-29 STOPS 


913.0 


590.0 


1.55 


-3.3 


GASKELL #6-22 STOPS 


1009.0 


626.0 


1.61 


+5.7 



Table 2 
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comparison with one generated for the same problem by another computer 

based method or the conjectured optimum routing system attainable. 

The ratio seems to be sensitive to the range and variability of 

the elements of the interstop travel time/distance matrix and the 

34 

grouping of stops as evidenced by the Thompson-31 stop problem and 

18 

the Gaskell #4-21 stop problem . Both of these problems had inter- 
stop travel time/distance matrices whose elements had a narrow range 
and low variability and the stops were arranged in groups. Moreover, 
both of these problems also had low ratios. 

This ratio appears to reflect the effect of using more than the 
absolute minimum number of routes required as evidenced by the Cochran- 
25 stop problem*®. The routing system developed for the problem con- 
tained two more than the absolute minimum number of routes required and 
the ratio was 3.88. 

The routing system developed by Newton's procedure for the Dantzig 

14 

and Ramser-12 stop problem traveled 4.8% more distance units than 
the conjectured optimal routing system whose total traveling distance 
was 290 distance units and the ratio was 2.10. However, the ratio of 
the total traveling distance for the conjectured optimal routing system 
divided by the lower bound on the total distance equals 2.0. 

Although the ratio is a relatively poor predictor of the degree of 
optimality attained by a routing system developed by any method, it can 
still be used as a criterion for terminating this computational pro- 
cedure. Knowledge of reasonable ratios which can be expected for a 
routing system developed for a school will be acquired after experience 
with this procedure has been gained. 

Although this computational procedure has proved to be efficient 



with respect to the quality of routing systems developed and with re- 
spect to the amount of computer time used to solve ten problems from 

6 

the literature, no problem involved more than 45 stops, the Boyer data . 

A thorough evaluation of a heuristic procedure requires that its 
efficiency be also examined for problems of the size which it would be 
expected to handle in the real world. 

Thus, 18 problems involving up to 120 bus stops were created from 
tables of random numbers to examine the computer time required to 
develop routing systems for problems involving various numbers of stops. 
The elements of the inters top travel time matrices ranged between 2 and 
20 minutes; the elements of the student load vector varied between 1 
and 9. These ranges of values were selected because they satisfied the 
characteristics associated with the school bus scheduling problem. 

All problems were constrained by a student riding time of 45 minutes 
and all problems were run until Algorithm A was exhausted, i.e. no 
further trial routes could be generated. A summary of the results ap- 
pears in Table 3. 

By observation, it appears that the lower bus capacity for. a pro- 
blem involving the same number of bus stops requires slightly less 
computer time. This is due to the fact that the individual routes 
contain fewer bus stops than the routes for buses of larger capacity 
and the amount of time required to improve each individual route is 
decreased. 

Moreover, it seems that the quasi-optimal routing system is 
usually selected from an iteration whose trial route was generated by 
Algorithm A. This evidence also appeared when routing systems were 
developed for ten prj^lems from the literature and thus bears out the 
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LY OF SAMPLE SCHOOL BUS SCHEDULING PROBLEMS 
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Table 



conjecture that better trial routes tend to partition into better sets 
of routes. 

Since, this computational procedure was able to develop routing 
systems for problems involving 120 stops in approximately seven minutes, 
it is considered to be efficient with respect to computer time usage 
and therefore would be of practical value to a school district. 

Finally, to examine the worth of Newton's procedure with respect 
to a real world situation, it was applied to four schools in the 
Williamsville Central School District, a rapidly growing suburban area 
in Western New York. 

This district maintains a fleet of 78 buses which service ten 
elementary schools, two middle schools and three secondary schools 
located within the 42 square mile area of the district and twenty 
private and special schools outside the district. Furthermore, it is 
anticipated that ten new schools will be added to the school system 
within the next five years. Therefore, the Williamsville Central 
School District expects to be continually faced with a complex school 
bus routing problem and expressed a great interest in using a practical 
computer based method for designing its school bus routes. 

The transportation director and the assistant to the superintendent 
of schools requested that Newton's procedure be applied to Academy 
Elementary School, Forest Elementary School, South Senior High School, 
and Dodge Elementary School. The first three schools service the 
established, densely populated part of Williamsville. Dodge Elementary 
School is located in the new, sparsely populated section of the area. 

For each school, the school administrators designated the bus 
stops on a large map of the area and assigned students to the stops 



from census tracts. Each bus stop was assigned two labels, a road code 
number and a map number . 

The interstop travel distance matrices were then developed manual- 
ly by using a map reader. Two people developed the matrices for the 
four schools in approximately forty hours. The elements of the inter- 
stop distance matrix were then multiplied by a time factor, specified 
by the school adminis tratuzs , to convert the distance matrix into an 
interstop travel time matrix. lhis time factor was large enough to 
include the time spent in servicing a stop. Since the school adminis- 
trators felt that variability in travel time was negligible when con- 
sidering all the routing systems collectively, no attempt was made to 
include the effects of this variability. 

The elements of the interstop travel time matrix for the Academy 
School varied between 0.5 and 14 minutes or 0.25 and 7 distance units; 
the elements of the student load vector ranged between 2 and 44. The 
large elements of the student load vector were due to the high density 
of the population in the area serviced by the Academy School and the 
lumping of points by the school administrators. The results of apply- 
ing this computational procedure to the Academy School data is 
summarized in Table 4. The routing system developed by this heuristic 
procedure requires 40% fewer distance units and 1 bus less than the 
current routing system used by the school. 

The elements of the interstop travel time matrix for the Forest 
School ranged between 0.5 and 14 minutes or 0.25 and 7 distance units; 
the elements of the student load vector varied between 1 and 30 
students. The Forest School is located in a high population density 
area. The results of applying Newton's method to the Forest School 
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data is summarized in Table 5. The routing system developed by this 
computational procedure requires 38% fewer distance units than the 
current routing system used by the school. However, both routing 
systems use 10 buses when the bus capacity is 67 students. Because of 
the lumping of bus stops, the routing systems required more than one 
extra bus over the absolute minimum number of buses required when the 
bus capacity was reduced to 57 and 52. 

The elements of the interstop travel time matrix for the South 
Senior High School varied between 0.5 and 22 minutes or 0.25 and 11 
distance units; the elements of the student load vector ranged between 
1 and 41 students. These large student load elements were due to lumping 
of bus stops by the school administrators. The results of applying 
this computational procedure to the South Senior High School data is 
summarized in Table 6. The routing system developed by Newton's 
procedure requires 25% fewer distance units and one less bus than the 
current routing system used by the school. 

The elements of the interstop travel time matrix for the Dodge 
School varied between 0.5 minute and 35 minutes or 0.25 and 17.5 
distance units; the elements of the student load vector ranged between 
1 and 27 students. All the routes for the Academy School started from 
the South Senior High School, all the routes for the Forest School 
started from the Academy School and all the routes for the South 
Senior High School started from the garage. However, the bus routes 
for the Dodge School can start at any of six origins. The results of 
applying Newton's procedure to the Dodge School data is summarized in 
Table 7. The routing system used by this heuristic procedure requires 
17% fewer distance units and two less buses than the current routing 
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Algorithm A was exhausted 
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Algorithm A was exhausted 
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Algorithm A was not exhausted 



system used by the school. Since this problem involved 96 stops, a 
study was made to examine the effect of limitations upon the number 
of trial routes generated. By observation, it can be noted in Table 
7 that increasing the number of iterations allowed for each case used 
considerably more computer time but did not usually reduce the total 
number of distance units required by the routing system. 

The application of Newton's computational procedure to four 
actual school bussing problems demonstrated its effectiveness. It 
was able to develop routing systems which were superior to the routing 
systems currently used by these schools with respect to the distance 
traveled and the number of buses used. Moreover, it was able to 
produce these high quality routing systems using an operationally 
acceptable amount of computer time. 



VII. SUMMARY AND RECOMMENDATIONS 

The model developed in this dissertation is a general one which 
utilizes variables that are applicable to all school systems. It 
represents a significant improvement over the manual and other computer 
based methods available in that it generates routing systems which are 
efficient with respect to the total mileage traveled and the number of 
routes required using a minimal amount of computer time. Moreover, it 
routes buses from school- to-school in addition to developing the in- 
dividual routes for a school. 

A review of the work accomplished during this study leads to the 
following recommendations for future efforts in the automatic design 
of school bus routes: 

1. The manual calculation of the elements of the interstop 
distance matrix is tedious, error prone, and time 
consuming. Therefore, a great need exists for a 
computer based procedure for developing an accurate, 
non-symmetric interstop distance matrix involving 50-120 
bus stops located in either densely or sparsely populated 
areas . 

2. A computer based procedure for determining a student 

load vector whose elements could be constrained by student 
walking distance and size of load at a stop would be 
useful. 

3. In order to avoid the development of a routing system 
involving a combination of full bus loads and half filled 
buses and a combination of maximum length and very short 
routes, a need exists for a method which would balance the 



sizes of the bus loads and the lengths of the individual 
routes. 

4. Since government aid is usually given to school districts 
for every bus load of a certain size, it would be useful 
to be able to specify both upper and lower limits on 
bus capacity. 

!>. School districts are interested in the cost associated with 
a transportation system. Therefore, it would be useful 
to be able to develop routing systems with respect to 
minimizing a cost function which would include the effects 
of the total mileage traveled and the size of the buses 
used . 

(>. A great need exists for u computer based procedure which 
would optimally assign all the routes which a particular 
bus would service between the time it left the garage 
and returned. 
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APPENDIX 



The Appendix contains a listing of Program BUS2 which generates 
bus routes and schedules for a multi-school system by means of the 
method described in this dissertation. Program BUS2 is written in 
FORTRAN IV for the CDC 6400 computer. 






■ ' : \ 

•-* ftj-m IT N OD ^DDO'OH(M^OF4fn^ll)^) 

o o,q' p o o o o H»^Hi\j(V(\jfv/ irtriotnihu) 

c o *oo oooo ooooooooooooooooooooo 

<\j m ^ it 
cc c c c 
a <ra <r or 

cc C C C — ‘ 











4 « 


4 # 




• 






CO 






Uf 






* 4 


« 4 




cr 




<— — 


» 






►- — 






* * 


4 « 


4 


UJ 




• — QC 


o 






Z c 






4 * 


♦ 4 


♦ 


H* 




a u 








C c *4* 
















U. ^ 9- 


<f f 






cr (vi oj 












4 




9 - uj a — 


• o 








— — 










— 




9- 5 ff 


a 






<r ~ a 


•X c 










. X • 




♦ X C • 


-» 






3 *-« 9- 


m cm 








4 


(o cr 




Hif o 








cc a a 










4 


— •“ o 




• w 


•» 






M 












in o 




or • •< 


— — 






X 


X CD 










if H .10 




o — — x*-< o 






a » *04 










^ Ul 




CD CD Jk 


♦ *■* 






<* — — 






Ul 






W w O 




to — ^ 


Ul > 






UJ C O 




9- 


UJ 




l-<n o 




DO J * 


< 






(T, fir 


«► » 






X 


♦ 


ax -J 




X U if — 


• CD 






a 






a 


»- 




2 X O 




o r> a lj •** z 






o s rr 


it in 










•-* 2 Z 




O J w • ^ 


4 






u =3 uj cr m 




vC VO 




X H 




m U — CD UJ ^ 






2 *- 


M W 




^ CD O 


H- 




• w 




a 2 • u» <r 


^ — 






o x 


H» H- 




X O < 






•» — 




J£ W < 






u. •-* 


a ir 




o a u 


cr 


« 


— X — u. 




a u. u 


• 9* 






9- * 


2 




<n oa 


2 




in n o »- 




^ ^ o • 


• O < 






• 


O • 




-j ? 


*— 




*vi » «->4 X 




cr —yr — 


— o 






(O O — 


^ — 






O 




w o ^ O 




hz ja x 


o a a 


2 




w(*l O 


in o 


• X -I 


X 




co-* a a 


— < 


M O • *-* 


w o ^ _J 


O 




> ~ <**► 


•> cr 


—4 


aUJ JO 


* 


O — UJ ft 




* or G — 9- 


ff J fit o 


►M 




-v 




W 


O <1 O O * 


x <n a -J 


4 


31 *Ci^O UiO 


»— 




uj 2 rv 


If! 9— (Vi 


H O 


O Uf 




X D D (/) f 




o o o 


-i * a x 


c 




TCI 1 


X 


— 


w **f 3 


C CO 




^ x cr ? r 


-J 


C tf < J • 


— 30 


o 




-1 9— 


— if 


0112 


Z> 




I 2 ZOOO 


#-4 o *«» t>4 


* <n <n 


M 




c •- cr 


-J 


i 


c ^ < 






— j •»* 


c 


a u a r < <r 


— X 2 


u 




i- •- 






_J X 


X IT 


* 


♦ ♦ • J> 


IWOC2 >- 


U f 4 Ul 


M 








m 


a 


O 2 


# 


— ^ ^«k (3 •— * 


OH HU .2 hh 


-ui 


9- 




CO — » 


* in 




* 1 X 






Hoiniotf 


> cr o ~ — o 


2 




a o — 


— r* 


C- 


— 20 


_J UJ 




r r- 2 




a a <o •‘if 


<Kffh 


UJ 




cmc 


w • 


»— « 


C Us CO 


X 




♦ ir»Mcr 


c uj a ► — 


c a — <o 


o 


? 


i~ — rv 


CM C 


*— • 


fV 9- -J 






— > — < 


c 


U <0 c — u • 


-J z — 


M 


o 


(/) K «— • 


*-* ^ 


r-l 




cr cr 




m < a 2 


u. 


9— <1 *9 — 


CL O 




J- 


Ct 


w 


c 




c u 




- c c c c 




<r co cr — u rr 


• 3 C »- 


c* 


** 


uj a: lj iijh arc <»a 


♦-I- 




CM 2 9- UJ 


<n 2 -J c — • 


— 9- 


o 


u 


X < 9- 


cr cr 




<■ c 


O X? 




^ cr - if 


U- *"■ c u < a o 


<J •> to 


►H 


2 


►- •-« a 


if c 




O 9- 


uj a 




* 2 c a 


9— 


e o c c ►- — 


< 


cr 






y 


» 


. tf c 


> X 




♦ - c 




»- T JO > 


X 


u< 


< 


c » 


*> 




2 U* 


c 




— O •‘<0 9- 


6 


a do: h h if 


<o ~ * 


a 


9- 


2 ^ ♦ 


•> 


iT 




cr o 




— <vj — ztoactou. o 


*-i <0 — O 




H-* 


< O — 


o 


(O 


— < 


2 




(VjHOili 




x * * * 


JDOtf) 


Ul 


a 


r ir 


^ IT 


• 


c v- u 


»- UJ 




— • — x u 


u 


U U T — — 


a w •- 


X 




2 — 


^ rr 


onzo 


> X 




• c^Hooccr < ^ + 


» >f (Vi o 


»-« 




~ X — 


-J ~ 


r-t 


— o 


-J k- 




M O 




— 


— < <o 


9— 


9 


013 


r-» X 




1-0 2 


o 


* 


f\J *^0 XX 


X X »- 'H -Q O S3 ac 






u. 


<f a 


<£ 


if o 


> IL 




h i r ii* uj 


UJ U : UJ J cc u — 


w Q- C . 


cr 




X < X 


> O X) 


*-H »— » 


z o 






XXX < X 9- z 


(M ^ Z 


o 


<V 


O m r-i 


•-* if 


-J Jttn 


M 




S 22 X 


X 


X X 9- <1 XX 


O — UJ 


u. a x 


v> 








02 


X 




9— X> 


D 


DOOWHC 


x •> a 


3^ (VI 


3 


UkJ 






cn uj 


CO *- 




<02 2 


2 2 9- ^ 


-J 4 — < 


-J K 


CD 


12 Z 


2 2 


2 


T J5! 


►- 'r 


a* 


2 2 2 X 






<X9* 


UJ » 




1 - c c 


oooo ~ 


2 2 


4 


C O O X II 


II 


ii ii n — — ► 


•> » 


(M CD in — ' 


T 


. f* 




•-i 


t-i O U . 1 UJ 




•—» ^ 




<< — 


— * o Ul 


< 


<i 


5- (/• U 


1 / If 


<0 


u a 


^ «1 




<f <0 cr 10 < 


a 


e c u: — ~ a 


o — — a 


C J ^ h 


X 


0 2 2 


2 2 


2 


z ox 


UJ 




2 2 2 >- 




I M w 


— ♦ H < Z <t 


C 


U. U U. 


u u 


U 


1 K II 


K c 




u u u < 




a u < 


— < cr a 


»- C C X 


— 


-i r x 


X X 


x 


7T 09- 


<r or 


$ 


XXXX9- 




J hi 


c-DUJ5<<rcr 


a 


u. »- *- 




H' ►- U, 


►- c 


•* 


►— h •— cr u* 




< arc 


a cv a a 


u a u c 


a 


CO Q C. 


o o o o > <x 


V)3 




O O Q < -J 




h w ^ J Z 2 1 


xa x a 










* * 


« « 












— a 












❖ * 












a 










4 * 


4 A 












a* 










* 4 


« ♦ 


* 












a 


o 






OOOO 


o 


ooooooouooooo 


o 




92 



POOR ORIGINAL COPY- 
AVAILABLE AT time FIL 



»h (vro ^ ifi h n 

NflDoopHftjfnmromfom^tfuntnininintrnn^KoD^o^ftim^in^Kffij 
UVl/VlP aT >C C a £> *C a£ a£ *C >C aC C «£ vC nT VC < i N N S h N h> fw> c , 
o o o o ooooo ooooocooooooo oooooooo o o o o 






■■/ 



c 

o 

►* 

cr 

Ui 

a 

© 



a 

3 

o 



r 



a 

o 

o 



«r 

3 

Z 

2 ►-* 



»■> 










Z 




w QC 


cr i-i 






o 




o 




UI H- 


UI X 






c 








< K 3 3< 


K Z © 






►4 




»- 




O o 


3 •« Z 






a 




< 


e 


►-* C 1 it 


© CD 






Ui 




1- 


Ui 


U» 2 H 1- 


a: . h • 






a 




3- 


cr 


< 2 3 MS O 


Ct 0 * 










O 


cr 


i- ‘ o :t •-* 


cr o . 






c 




a 


UJ 


z z x i ct: 


3 »-» M 






z 




cr 


u 


lu H Ct »- IjlI »»• 


X ... b" 






M 




z 


o 


o c •— z -E UI 


Z < 






a r 




<r 


ct 


MHUJ^HN X 


• X/- V ’v ■ •-< 


Ui 




3 




X 


a 


2 o X J-- X 


O O Ui * 


-J 




O 








JO < | CL >* 


u ..... J 












cr 


C ►“ h- Ini HO 


a: o: m 


< 




e 




cr 


-i 


o cr c z uj x 


cr :■-••- :0 <X ~ 


— J 




Ui 




o 


o 


x ui ui ►-* *- ►- cr 


;Z; v;- «J > 


1-4 




u 




u. 


o 


ucx cr < *-< 


►*. a h < 


< 










X 


cr 3 o a h it 


i. u © ■ O <1 CD 


> 




> 






o 


-J UJ oox 


■■■>-,. U. > 2 


< 




a 




< 


cr 


IT or CJ H UJ Ixt »-i 


oc < 






UJ 




aa 




(vo2Zhi>»a 


v j-, C . W :i ♦ 


cr 




IT 




i 


u 


f-li H Jfi Or cr H- 


ui cr a 


Ui 










o 


M ||ih> < 


vV;U-;. X U >~ 


cr 




cr 








* rr » u c > u: x 


-j ^ cr ^ 


— i 


— * 


-j 




o 


cr 


— c cr c uj z x 


>.«CT£ 'L,. ■";;«■ i r cr 


a? 


•— 


c 




z 


u 


c\ o. r »- x u 1 


• K cc © 






c X 


< 


CD 


cj c 2 cr >~ x 


cr a a ac 


cr 


> 


X u 




X 


• Ui h c Ut UJ cr ►-» 


^W O,^UlU;; ►< Ui 


< 


u cr 




2. 


«-*"> (/)! hh 1— 


tV“4;v K O 


T 


cr 


IT 


z cr 


Z 


IV 2 3 U-f 


:. a.> r^. cj ,:. • x 


z 






3 




C IT 7 N O X UJ 


■'-■*. ^ M ^ X ^ 


3 




U. 




OD X 


cr 


X cr 3 O Ui b- X 


U. O X UI NO z 


X ♦ 


O 




Z X O 


o 


U OOOUOX 'r- 


.re- © ;.*■*.< a>*-r <r 




' © ... 




o 


c cr 


(L 


IT X ? 2 U QC 


- . o ... z X o 


-J 


Z X 


ac 


o 


z z o 




Z U IlQC o CO 


X O • U D P «- 


<t 


.., Z> ■ 


UJ O X o 


X 


cr so o o cr cr u. z 



,,£C : v , 

X ♦ ■ :•: 

3 l/T ~ II 

Z* ,-c *. '.' 



U .f ^ ',:< 

C 3T 3 7 

cr < X; cr <r,<,sr 
c nee a u c 
•Z S Il J Z S I^ 



^ It ? 



- X.MW » » II U.I 

<r o *■?•.« .x a: -< *-< a +- 

n > »- o ? 3 •!- — z ^ i 

. 1 / 1 w H. W2 tfl< « n ->-UJUJ3 * 

-rtf U‘ Lilt- » n 3 , X,: «► 3 3 3 O 

*- 2 II w ►-« 3 •- * Z X c C U. 

x x ; th a 
o 3 ir 2 k x c < a 

C/5 cr <TZCJ<>UJ 
a a c: a c cr uj _? c 
z z a z uziruo 



V: I 



(T yrjh* K II II < 

3< Z Z C: 

caa co •— i*-d 

O Z -J U O ^ X 3 



la. Z Z <1 
X «X^UJUJl 
It Ui UI Ui Ui O X 
iT O X X O 
X C U r; 

3£ +- Z L O 

Z- II u 

It II It 

it in oi 

Vm/> cr in in 

« y: 



K ~ X 
cr o 
' lx u cr 



X i— n* j z X X X 



k a 



CJ 



u u 



u 



u 



u u uuoooou 



93 



ERIC 



POOR ORIGINAL COPY - BEST 
AVAILABLE AT TIME FHM** 



O' o oi^rvi fn <* in vO n ©O' oiNMm^©©^© o o 
h-crcoa occcocococ oc ce cr a o O'O' o oo 00 s ^ o 

c OOOOOOOO OOOOOOOOOOO OO in 



COOOOOOOO 



© 

UJ 



c 

DC 

u 

o 



UJ 

0 



a 

o 

u. 

u 

WO 

a 2 
3 
uj o 

3 : 00 
2 : 

>- or 

CO Uj 

v* O 

O _l 
2 

UJ 
CO > 
W o 
AC 
X < 
►— 

a: uj 
1 — tr 
< < 
X ►- 
2 
u u 
3 - e> 
*-• a: 

i- UJ 

c 

UJ 

X UJ 
+- _J 

cr 
to < 

2 3 P 

< o 

UJ _J 
X -J 



o 

UJ 

o 



CO 

UJ 

♦- *c 
3 

o •* 

x •-« 



<r 

M I 

cr 



<0 

10 

<x 

a 

vr 



cr 

UJ 

a. 

* 



co 

cr 

x 



u. 

c — 

X 

& c 
uj <n 

CL -J 

2 ^ •> 

3 ~ 2 



<0 

¥ 



*-* 

o 

or 

o 

X 

O 



II 

II 

C 

CO UJ 
cr a 



3 (V 
X O 

1-4 o 

X I-* 
<t 

X • 

LO 

H 

to 

tr c 
< <r 
a u- 



< 

\C o 
O 

* o 

X ~4 
sC 

* 

IT ♦ 

2 

K W 

< U 

2 C 
fX * 

c u a 
u. x 2 



cr uj 3 2 ^ 

O 3 > D 
2 2 » •> 

3 0*^^ 

X 2 ►— O 

*-• M H 

X li || 



a 

o 

x 



> 

►— 

a 

w 

< o 

X UJ 
3 C 

X o 

O _J 



X 




2 


a -j 


u» 




O 


u. < 






a 




-J 




u. 


own 


UJ 






UJ H O 


> 






►- h- 


< 




U‘ 


U J (0 


ac 




> 


2 UJ UJ 


►- 




< 


-J <0 






X 


UJ <3 


0 




h- 


X COX GO 


*- 






++ *- 






O 


X UJ >- 


0 


2 


1 - 


»- ADO 2 


u? 






< Z < 


ac 




O 


2 -J 


►— 




U< 3 


JVC 


3 




X 


U. f M©Z 


0 


H 


M a 


I?lO< 


UJ 




3 O 


►-I < 


cr 


3 


. C H- 


H- ^ O 3 






UJ <0 


UJ 0 


to 


•> 


a 


O U 3 X C 


UJ 




in 


»-• 0 X 


►- 


3 


if 3 


0 ax c 


3 3 


» 


UI 00 


*- 1 - 2 CO 


2 


•— 




U W U' 


^-r 




3 c 


2 X cr 2 X 


2 a 


X 


2 h 


X X 3 x c 


0 


O 


H 4 


>- O* >- < 3 


u. 


2 


2 


co ^ co a co 


O CO 




H 


UJ 03 


X 


X w 


u 


cr »ticru? 


or ©0 


O 


c a 


M 0 H 



♦ o 

O 



n ^ m 



1 < ~ 

2 X 3 

» 

It X •-» 
< 

^ii?a 

2 — 1 c 
202 



< tr •> 
cu ru »-* 

X 

c c o 
0 0 2 



UJ UJ 
3 3 
2 Z 



IV 

(VC 



c c 
u o 

r ^ 
<v cv 



inuj c 

32 
LT ^ 2 
PU Hfll 
C K <x 
C?UJ 
o U- c -J 
Q X O u 

If 
CV 



0 2 2 CO >*2 2 

X *- 2 U. 2 _l X 

UJ CO » » •* — i _J lL •“ 

CD «- 4^-4 XO < I <H (\J 

2 to ►-* mx uj 

D D IV II X (\JXUJtt 

2 x o »- a 1 

wo <f Ul to to 1 - •-» 

II IILJIiJ 2 V) 23 M 

3 3 CO O C 2 

~ a k ~ z z 1 x -i a ui c 

3 (VWOhh Z UJ ^ 

* »I-Kh ^CC 3 

*-» ^ 2 2 tO «*<:*»- 

w CC^C Cli lLltliCli'C 
2 O O 2 U O H* a XX CO o 



N cr 
cv cv 



a 



00 u 



uuu 



uuu 



uuuu 



Id 



ERIC 



94 



kflitiMUTTihTiMJ 



POOR ORIGINAL COP' 
AVAILABLE at time I 



^nqo «*© (vi cdac ncd ^oh i\i ^okoo o^o «-«hm n hum 

*-*(vjm tv cvj tv tv tv cvt*>fo ro mnrom cn ro r» ^ «* * >^ * <* * m 

■4 • ■ ai mm* mm* mm* mm* m^ m ^ mm t _* . . ^ . . ^ ^ 




It 

3 



in uj ui 



x 

n 



a 

c 

x 



X 

H 

X 

H- 

< 

X 



o 

tr 

u 



o 

-I 



z 

ui 

o 

3 

H- 

(/) 



X 

►— 

a 

h- 

< 

X 

u 

o 



< 

X 

o 

cd 

< 



II 

3 

» 

t/> * 

tn -*•- 

< 

Gl X 

it: • 

W-* 

• II 

a: *h 

Ui 

a • 

* -»* 

3-1 

• • a- > 



tr -i 
x o 
<o 

X X 

u 

» <n 





c 




M O • 


r*. 


c in 


~ < 


tr 


UI 




I* w 


OC 


r 


X X o 


h* 


H 




tr 


a ui » 


o 


Ui 




o 




tL || 


U. 1- z 


tn * 


~ + OC 


X 


Ui 




o in 


3 


►-* tn 


v.^% 1-4 


o 


mJk 


• «*■* 


K O w 


UJ O -M 


o r 


^ 3 


Ui 




> tn r 


t rr e 


z e 


M 1 


<n 


X' 


» 


WS QC 


a • 


» * UJ 


Of 


x tr 




0C - 


1- > 


C 2 


— •#-. cr 


h- 




3 Kw 


Z • 


it 


Ul 


M 


I-# 


X — 


CD <t ^4 


UI • 


z w# <n 


x > •• • 


X ^ 


H 


Q Q 

or ui < 


2- .. - • ' 

»— I’U. mm* 


a -* 
< 


•: uj 

r- • tn 


>* 


H 




c z © 


a a *- © 


i- >n 


II 3 


t ft 


tr 3 




u e -j 


•w a o 


» 


3 *w a 




h © 




w 


a x 


2 rs 




u 


S? O 


<r 


— tn » 


— <-2-. 


c n 


• c u 


o 


X 

1. 1 


c” 


UJ V) 


Ui iw 

1 A ■••■"- 


M 


^ < o 





u 




u tn 


1 » 


<x 


► 


1 




r r 


►- 


. -y:.C 




s 




►4 < 


■ " w> • 


<■ 


X 


X "■ 




H u 
Ui 


■ X 


a © 

UJ mJ 


w 


cr pw ; : 
ui x 




c tr 

M 




h~ 

o •— 




CL 2 


2 


OC 


z « 


<z 




CL 




*- z 




X ui 


o 


- X 9 


■:♦ »— 


u 




< c 


o 


■ . 




X 1 


AIO 


X 3 


© 


UI 


, ■: 3 


X 




O f--' 


— • 


X H 


II « 


>- ru n — • 


in 



UiUJ 

r 3 



tr 

10 M > 



cc c •» it 

CT C « h* _j.l~ 

< < J 3 IT CT Z 
I C X U> 

0 -10 < i: 
'3 J h \ X Uv 
tt <T vJ 

a ~ az 

X X uJ 

> 3 3 Z ©> * 
XX © O V 1-4 

H M.M © M 
X X 0 ^ ? II 
< < IT •- 



:'.*■* 2 . < :■*•**• r 

M tr sn- 4«*3 
C H ll; ••-.> 

<r X ,2 It •-« 

ur © u c : e ^ 
tr u p o r o x 

c r 4 

^ * St 



(VZ 5 
IT OC 



in ui it 



a© - 2 - 

»W* •■* © UI •— < 

H H K C C C l-~ 

x z < < <7 or 

© e c u* u c: u ; c <s 

oo o cr a: a a : O -J 



II • X X < 

in 

— ii n 



a 

* 



.e 

IT 



c c x a 

<< XX 

C U < <3 

-J tr x x 

tv 

IT 



» ^ H' 

IT 2 

■*-r.. 

w 

uj in 
c a 
< o 

U I- C' 

QC WO 



0*0 

CM> 

O'O' 

0*0* 



ii n ii 



© 

o 

T 
U 
<A 

LOT 1 

o o 

<>Uv 

o <x 

'(W» 

ui n < 

3 - K.C 

? — 

w 2 UJ 



XX X 



eva 

0X3: 

* 

IT 



• tr -j a 
— c u 
*-< z ^ <r 
^ » x 

X 3 

c ir > z 

0 *-!■;'•: .;;• 

Z X 
• *- <► 3 
X 

•> 3IM H-I 
-*-• ► w 2 
^ MK M 

' - • Z 

^ X ~ 

wow n 

x tr 
ui ^ < 3 

k -j cr 

w ■ . 2 

ar - ^ 

3>,' X 

: .^CV 



in ° 

S £ 
00 *< 
# 

8 W 



2 < 
o 

2 “! 
O » 

o < 
2 | 



U 



o o 



oo ooo 



o 



o 

ERIC 



95 



tn in in in in in in * 



*- (v n in 

r\jn^ininini/nnmi)NCD(> onivjn «*> in <n 

sC \0 nC vCNC'n'CNC'O^X^N^N^' h- I** K 



•-• cm cn in 

»H H ^ H ^ 

%c. wc ^ >n 



* * 
* « 









♦ * 


* * * 






















in 




















X 




















— > 










*— t 








* 




















♦ 


a 
















l/> 




z 










_J 






z 


Ul 


H 






m 




o 






o 


z 


1- o 






UJ 




o 






M 


a 


z z 






►-» 




X 






H- 


♦ 


UJ M 






-J 




u 






a 


2 * 


c u 










if) 






HH 


< 


z o 






H* 










X 


X 


uj u* u © 






Z 




> 






1- 


►— 


a in c a 






UJ 




flD 






m 




UJ 3 Ui 






X 










UJ 


X * 


o cn ni 






UJ 




c 






a 


U * 


z in in 






-J 




UJ 








►— 


H4 UJ UJ c 






UJ 










tn*- 


< 


H- U 










o 






►-X 


UJ 


l/>3 O UJ 






«— 




-J 






Z CD 


X 


•-> O X > 






in 










Hr H 


<5 « 


x a < 






UJ 










< h- 


♦ 


X x o 






-j 










a 


in x 


OlLUi X 






-j 




if) 




X 


H-» 


3 U-* 


hOIKUJ 






<r 




UJ 




O' 


inx 


OP ID 


X HON 






X 




in 




a 


Z UJ 


be U* 


3 X Z 






<n 




3 




a 


o o 


H- * 


uj cr < 










CD 




a* 


o z 


4 Z * 


a x z in 






X 








X 


3 


H 


z z h* uj cn 




HI 


o 




U 




• 


e 


IT 


33 a c < 






H-4 




o 




o 


Z Z 3 Z 


O Z 3 Q X 




r- 3 


X 










•-* o 


X < 


CO C 




o 


3 




cr 




4 


Hh2 * 


UJ 




♦ X 






u- 






1A H- 


h- v # 


or T C 3 3 




•— 


Z 




00 






33 X +* 


UJ h- UJ 




r- z 


H 




J 






i 




3 X Z Z 










3 






o c 


n 


*- C C U* X X 




» 


z 


•J 


z 




a 


< IT 


IT 


-? z ♦“ r r 




C H~ 


X 


c 










~ 3 * 


4 < J J J 




Z 


3 


o 


X 


3 


<C 


Z < 


H* AC * 


u < o o 




UJ 


_J 


z 


3 


w 


X 


UJ 




in x z oo 




** X 


a 




X 


•- C 


c 


X z 


w 


r hh(t 




u 


o 


UJ 


M 


X «r 




•?• *-< Irt 


x c h* in in 




X Z 3 




o 


X 


z o 


N* 


< Z> IU 


ZU.HI z z 




Z Z h UJ 


u 


r 


<T 


-1 




UJ *— 


X x * 


~t O X <x < 




X 


o 




X 


» 


-J X 


e d 


X liv # 


x c -* u? u> z 


» r H 




c- 




CM 4 


3 3 


z z o 


< X 


Z 3 x X 




CM CM t in 


X 


z 


II 




x x in 


H-# < 


x it 


II O UJ » 




UJ 


^ UJ 


< 




« X 


D< J 7 r o 




W X HI 




II H- fil — • — J 3 X 






3 


<n x < 


<U K 




— H- h- ^ O r-4 


o 


f- 3 -J 


•* X 


in 


*— r 


3 in 


-j 




* 


—» X X 


II 


H H3 •* < 


*-♦3 


CD 


r 


© JUJ 


H — 


ii -j a 


# 


3C O U It II 


II UJ 


N »-» X 


W Z /UJ J 



~ 3 *1 

V* H r 19 Z -J II 

3 >C MT CD 
<r ? 7 k r x 2 

x r> 3 

c me it 

r jq -J 



3 uj 
if u . C, 
DUiX 

a me 

7 (/)< M 2 3 
C J i i »- ^ 

U<<< 7 <h 



w a C 
V 

3 ii 

a: 

$ x cn 
• * <3 a 
x -j 



X — — * 
H-. 3 3 

or ~ w 

UJUI 

<x o c. 

7 



3 3 o 

^ 2 

in 3 •-* « 

* ~ 

UJ 

c: c 

O X 



z in 
e cr 

u -I 



K U 

c •- 
O X 



rw ii u 



H IIZ U; 



z 

o u ►• 

o ►« X 



c 



• * * * 

* * * $ 

* * * * $ 

$ $ * * 4 

o o oo o 



IT 

*C 



uuu uu 



Z O O Z 
HU uc 
x z z u 

a ^* 

* h- 



u 



o 

a 

3 



o 

ERIC 



96 



»0h-a5 0 'OHftifn h ®o» o »h w m ^ in ^ n ® ooH(vto^mvfiN m o 

Tf- t*- OC QC QC OC CCOC QC 00 00 00(7* (7 (7 C7 O O' (7 (7(7(7CCCCOOCCCC 



2 
1 U 
U i 





H l-» 




UJ 












CD 








2 2 

M H 




© 








© © 




2 




















a oc 




-J 








o o 


.. x ■ 


UJ 












> 








f- ►— 


X 


.< 








<t < 


M 


QC 










tz 


t- 








►- 1- 


1- 










tr oc 




>* UJ 








< < 


X 


OQ >— 








1- K O 




3 








(/> tr ui 


Uu 


O O 








cn 


O 


UJ QC 








c/> tno 




J- 






h- 


UJ. O tiJ 




3 >~ 






2 


H Ul KUJ 


H 


X 2 






UJ -.; 


r if ra 




< 




2 2 


x :... . 


C 3 o <r 


3T 


QC 




»-4 M 


UJ 


QC X 


O 


7-tU 


o 


^ X X 




uj tn 


QC 


2 O 


00 




UJ 


-J (X — J UJ 




C 




2 * 




-J < JH* 


C 


u a 




*-* >- 


X 


<£ <3 


1— 


. c 


o' 


x ^ or 


3 


tr o 




o t- 


GO 




X 


2 u 2 a 




2 tr 




♦ V 


M 


UJ fr- UJ 


QC 


3 




w X 


2 


i r 


O 


C 1- 


IT 


tr tr 



z 
•— « 
X. 

H ♦ 

~ (T 
-J CD 

© -J 

c 

2 II 

UJ IT 

car 

^ -l 



2 

U; ^ 

X X 

I 

>- z 

Xt 

♦ 3/ 

rv * 

W M 

tl ^ 

* x 
o 3 ■ 
tr u uj 

— s. 

:Uj'rr — 2. 
0>-5 hi 
3 ; ■ *1-1 

e: ' >-.z 

U.r O ^ C 

q: a x u 



3 0 3^ 1 OC tr. CO 

a y q 

c c w .. >c cr *-* 

2 V5 Z tT ( w tu u» 

D O DDaa 3* 2 

c or c t c k c — »- 

aczff y z < — -j uj xx 

•H < X Z O X z 

QC; ;X *;X ; • • * U. >- I 

UJ U . *-r U. * •""» O * 

3 0 3C O O *■« ^ O (Vi 

0 2 0 2 11 X 22 

-j < -J «x n o o <n 11 3 

* X *r* N * 



UJ 

3..~ - 



tl 



•-* «v 

20 O 

o or cr 

U-J -J 



c 

Q 



3: II 

O ■;*■ 

OC 

Ct 3 
U* * 
> —t 

c c ^ 

X O 2 



Ul QC ►-* N- 

3 o h- 

2 0 . 2 0 * 

M M II ® 

K CC 
??+-»2 
c a i- cr 

uik 0X0 



UJ o 
2 U 
II 2 2 

II i-i 



- X 



r if 



0 2 
U H~ 
2 X 

IT 

cr 



2 W 
c c 
u *: 

a 

x 



3 3 X 

XT'*-' 2 

?xa 
>-*<H X 
X X < 

r h i 
w ^ X 

Ul» Z *-» 

» ♦ ' ♦ O 3 

cn in cvi — « 

X X -* 

-J -J II X 
3P 

II II 0 3 II 

QC UJ 

^ ^ —' 3 

*- — iv e. Tz- 
W w O O* *M 
— (V X • ♦— 

O OO •’ 2 

a: a u: tr v.c 
j- J oc p x u 

c 

a 



r> 



i*> 

o 

2 * 

rn 

- •»'.■ > 

•<Vi 

n 

— i 
3 31 
— I 
I Uj 

m c 

h 



r- 



c u 
q : 



r 

a 



o ouo 



o o 



o 

ERIC 



97 



If < m(K*J) - min 
97 MTN s M(K,J) 

loi Continue 



o^*ruf»i«>-ifivGL^cc o»o no cv cm m in sc ao hm a> 

ojoj dj cv rv cv ou tn m ir in in in in in »c *£*£>£> c »o >c so 

(\j rvMM (mm rv oj (\Hii m <vj cj r\j(\j (\j fu t\i m oj a mcm mvi (v<\MVi(\i cvjm m m mm 



<c 



r- 

a> 















c 


















UJ 


















X 














UJ 




o 












— > 


2 




1- 














o 




<n 












2 
















• 


U 2 




> 












r- * 


C U 4- 
















(X II 


X X 




in 












o 3 


H- 2 4 




3 












<1 :z 


m 4- 




O 












« 


o in 




M 












o 


O H 




> 












<n m 


*: in 




iU 












H 2 


X UJ 




X 












C II * 


» jk 




c a 












FM 


X 03 




UJ 












x*-* ii 


^ o o 




4- < 












O #■* 


2 x a 




< 1- 












» 


1-4 U 




-J <f 












UJ * 


X U)J 




3 O 


in 






IT 




a h ^ 






O 


h- 






r^ 




< ~ 3 


O X < 




-J -J 


2 










H<V • 


u. c 


M 


< o 


UJ 






#- 




o *-• 


4- UO 


p-4 


o o 


X 






2 




2 00 - 


O 2 




X 


UJ 






Ui 




c jr 


LU <X 


w 


2 U 


a 




z 


X 






X M 


in 


ui in 


*-• 




J— « 


ut 




*■* 


X 3 


3 


LU 


3 




— X 


►» 




*«♦ 


UJ Z UJ 


X 


X o 


O 






< 




* * — 


M 1/) 


X 


2 


UI 




X * 


h 






X C 3 


<r 


U fc— 


X 




fr-4 


<n 




_j — *• 


O M 


X 


> o 






X ~ 






o — 


X UJ 




< < 


in 




H 


u. 




O X 


o o a 


\ 


X UJ 


3 




« * 


c 




X C 2 


2 • < 




X 


Ct 










0 2* 


3 h- 




» 




in ir 


it 




in — 


c < m 




<J X 


o 






3 




» II 


CD UJ 


o 


l~ c 


2 




c a or 


f- 




Q: p- 


cr or 


■© 


u u 


<r 




o z x 


< 




C 


fUX 2 3 


•-4 








#— ♦— <r 


4- 




U If * 


— ' U ‘ X 




in 


> 




Z 2 XX 


W 






3 4- 4- 




c < 


1- 




»»—•—» 






C N* H- 


♦ CG^ 




2 


4— 


4 - 


a: s v ^ 


o 




l_ ^ 


f*- U c F- 2 




3 


-J 




o 


*- 




< * LU 


*- b- X ~ 


4-4 


o o 


1-4 


in 


<-» ♦ 4 ♦ ♦ 


w 




a — o 


ui(/)^ocr 


>— 


cd in 


CD 


3 






4- X 


* C IT 1 - c 


cv 


*— 


< 


X X 




c 




c ~ 


^ < U? 3 X 2 


c 


x c 


_l 


U X 


in « M M M 


XL 




2 


f*h X 1— 3 0 


3 


Ui 


1-4 


in < 


O'—'—* W» ^ 


o 




30 


'jJ 3 X •* 


-J 


i X 


< 


2 X 


^ fv h CV 


a 




C X • «-• 


•-» > C < 4~ #— 




c o 


> 


o o o o 


O If 




x -j cn 


<a x < 


W 


_J 


< 


» 44 


— - OC CD CD CD 


»— 




x if 


ii 




UJ 




pH 


-1 -1 J J 


o ~ 




X ^ -1 


it 


If/ 


a 


in 


X 




UJ 




UJ X 


1 M 




~j < 


3 


n 3 


o n f*> w it 


> -J 




3 ~ — o 








X- 


in 


1 H Ui 


O 111 O W <M Z 


X — » 


^UJN< (VI 




OHM 


-*•* *■*> M ^ 


UJ o 


-4 


-J ^ — 


o*r m 


*: 3 


UJ 


ui 




T } ^ ^ 2 


a x 


X 


1 


4C • (V 


'!► 2 


ir a c 


ct 


h >n u 


^ ^ O ^ ^ M C 


m ut UJ UJ 


4-4 f-4 


4-4 M 


OC\J < 2 


< 


ru 


w •— «v i- •- cv »- 


K ^ 


b- 


KHH K 


w w 


H- 


H P- d ►- Q. 


X p- X 


o c cox 


in uj 


X 


H H H (/I 


c 


<;f 2 


UJ 3 ? 


r 


3 3 


iitraca o c 


u_ c 


c 


ctitr u 


u. fie C' 


IV cz 


OOf U‘ 


c 


trc it 




joo<a x o »-• o 


it o rr 








<C *F- 




IT 






C C r- 




*— < 




<v IV 


(V 


<V S 






H H M 








»4 f -4 


P -4 


f «4 








o 




u o 


ooo o 




o o 


u 






98 



(V »*« M 

O'a'O^ooooo^w^iDsos cd o o <h m n ^in^sco^OHivm ^ ia vd n 
f^f^f^h-Kh-r^h-acccccccccoDccaccflCO'a-cr o & cr i> a 
A*rvuv)WMw»vjw(\)(\i<vjrwA)MMAjrurvn\ifu(\»fu(vK\iN(vj(\jfu(utvj(\j(VMVfvj 



CVJ 





c 


X _> 




z 


0 0 




<■ 


0 

X X 


• 




U 


0 




in 


Ui 




z 


0 : 


z 


l-H X 




•— 1 


IT O 


3 


0 


•-* u. 


O 


1-4 


QC 


Ui 


~ac 


0 in 


a 


0 


UI 




ui 


Ui *~ 


in 


c X 


X 3 


Ui 


c 0 


1 - O 


ir. 


¥ X 


3 a 


3 


► u~ 


X 


CD 


in 


-JO u 


X 




0 0 


«$■ 


OD X O Z 




z <r 


x >-» a 




•» HP 


0 cn uj 


» 




sf 


OC 


a x 


H 


jitroD 



X — 

x in 



x »uo i z 

X !/)•-» oc 
►” rnrtL> 
HCC J 3 Q.V) 



in 

m 



in 

/n 



oc 

OU Ui 


<n * u; 


a a x z 3 u 
mo < </} h <ni~ 




•> 


— c 


cr a _ 


- C 5 U X <5 




IT 


<r . 


Z O GC 


-J ^ in a ui uj 




n 


* h 


cr <* 


* < >3 U ’3 OC 






a a 
<v c 


• or j 
7 u *- 


<r <cLhhe 
r. 2 a 






** X 
in 


3 x <? 
in tr > 


^OZIlCUUJ 
** * C h- O X 


X 


X 


*** 


*-4 < 


♦ H* X E- 


0 


X < 


<r 


^ ^ * (T 


e < 0 a .j a 


V* 


e 5 


r -3 ~ 


m k in w u* ui m 


z 


z 


3® 


X Ui 


~*cr z x > 3 ui 


k 


m 3 


O A i) 


3 * E < X > 


« 


•> «•*» 



Z 

M 

o 

M 

a: 

o 

> 

oc 

a 

tu 

u 

H 

> 

a 

UJ 

CT 

Ui 

a 

< x 
o 

IT IT 
-I Z 
O 



G 

Ui 

in 

> in 
-J; Ui 

-a o 
<z o 
3 a 

O Ql 

> cr 

— i 

Cl o 

2: o 

X 

o 

-j in 

o 

c u 
a: o 
1; 

in a 
u 

X CD- 
U X 
< 3 
l*i Z 

0 x 





X 1 O O* O 3 


— a a 3 a 3 a 


*- -*« ^ 3 


cr 


X — 


«- \ 
►"4 U 




z 0 a* 0 * x 


^ in 0 z ♦- z u i 


3 •> 


z 


u *-* 


0 ) 






z z tn 


II II * »-♦ 




in h 


O') X 




» X X tf> 


M < II || 


1-4 ^ 


11 


11 


Lrul UI 




3 » » » 0*4 


1 a 


3 in 




-j 3 


O »— 




X *- v© >C «C 


a ►- z ~ ~ 


m 3 




m .j ^ 


O Z 


111 3 Ui •> 


0 a <3 3 


3 £ h y 3 til 


^ < 3 


UJ OC 3 


3 


tr x www ^ 


Hzconr • 


Ncirat 2 3 


nr 


*H F- —*■ 


3 CL O 


2 


•— h •— 


2 h hi- a 


ir 2 it z 


a 2 


a ^ r*. 


2 U O 


»-4 


SUiUlliK 


O 1 3 ~ Ui 


<-» 11 *-* 11 »-» 


Ui 1-4 


0 rn <-4 uj 


M h* 


K 


''C KKK 5 * 


> H V y J (/* £L 


■ 3 h- 


a ♦- 


a 


Kft II II 


Z 


T M M Ml X 


m uj jo 3 


X X C 7 D 2 


3 


Z < 


c 


u a cr a cr c ur 


0 u u x < a ir 


cr <s c u <* a 0 


If c 


C K c: IT 


cr i- *•—?»— 


0 


MUJ 3 2 2 U 

1 fh 


OHnwu 2 z 


0 x a •-* z z 0 


zoo<ozocnxx 


< 


f- c 


x 0 


cr ir 


r^ 


0* 


;•- r 


A. 


A O 


c\: a 


rr r* 


rr 


fr 






H O 

IV 


r -4 H 


•H «-4 


^-4 


p -4 


f 4 ^4 




0 


0 0 u 00 u 






O 


0 0 




99 



o 

ERIC 



tVi ^ ^ N CG <V s|- \C CC (V r ^ Ct (V ^ CC (\j ^ sC 

~ I X«'X I X'X < £~ ® “:' ® J?® J»ooqoohh-,^n(»i Mt- in >C CD tb o 
O^ff^^CTO'O' O'^ <T 0 k <^^C 0 'OOCCCOCOCCCCOCCCC-i 

cvcvcvcvfiicvrvrvojfufvrvrijcurvj wrnnnww nronn^rom 



if 

"5 



«x 

c 

o 

o 

fr- 

et 

Uj 

a 

u. 

c 

cv 

H* IT 

0 — 

c 

1 » 

IT 

frH H 

X PO 
Q • O 

in © 

Ui ^ -i 

h-- ^ 

<r » 

c 

Ui cn — 
5: * 
rr uj 

u-'K 

H 

z u a 



<v 

IT 



X 

CE 

I- 

< 

X 



X 

> sn 
o cv 
o" © •> 

^ »H © 

X 



X SO 





X 




X 


X 


o 




o — 


o 






in 


in 


z 




z cr 


. z 






o 




fr 




• z 


* 


w+ 4 






H 






fr 




n 




U 


H 










-> 




n n 


"> 














fr M 


fr 














-> • 




•> 




Gt ^ ^ 


* 


fr-. 




u cr *— 


w in 






a uj > ^ 


w 


V 




r a < > 


< X 


O 




tr o <r <* 


cr z 


cr 




Z tr z cr 


-i X 


z 




X Z X z 


s* 








H ^ 


fr* 






— — ^ fr or 


j-, ^ 




+> ^ nr j-». ri. 


X o 






X X 


©-» © jy z 


© — * 


-fr i/I cv Jl -*cn 


(X (v cr 


cv cv 


m 


cv rr cv rr 


© © © fr fr 


© © 


© 


fr © © •> © © 


iHH O ^ 






© *—»■—»©»-* 


X 






X X 


* • n 


fr fr 




H fr fr f — 1 fr fr 


sC sC 


so vC 


c 


vC nO C 


UJ ^ fr-« 


UJ 




W w 


" — O ~ 


W W ^ 




W W W 



T 

o 

l/> 



</> * 

in ~ a 
< o 
a. z z 

^ • i- 

H Q 1 fr 

• IjXH 
QC •-• < 

Ui X II 

Q. fr 

* ~ • Mr 

r» ^ » Qj • 
X (/)wX« 
ir^ <h 
J^XX- 
o cu 

0*^*0 
x in ^ ec 

os<: ^ _i 
in z 

» » * ^ 

az^n 



iu uj < 

V- K 2 

^ *— ar 
cac 

23 U. 

a 

IV' 

© 



f 1 2 . H ^ Z K 

in UJ Ui H Ui <* in UJ LU fr-r UJ C 
fr-H“ ►- 2 K K ^ ^ X 

h m ? mq: ^ m* jr 

caccacccra 

0 3 3 0 



»- > y i- 1 - 



r* 



3 U. 0 3 3 
'C 

rr 



cacctrcaa 

02U.33U.23 



<* r- 

it rr 

© 



Cl 

rr 





c n » 




U. ♦ z *-• ++ 




— » fr 


fr 


fr cr 


*p» 


< vC II — O 


“D 


— fr n z 


CO ^ 


O^ w 


O Kf 


fr (T II fr C fr 


X o 


O' MM ^ <H 


X X 


sf C O 


< X 


^ •* a if 


5 «? 


a : — o 


X X 


fr © M Z M 




O w • 




Si UJX IH fr fr 

fr- a c n — 


fr «*•— 


<r m *-< h 


X 


X ^ H* «l Z ^ 


J> H 


O ^ • fr M 


cv 


m x x ■«* •-» o 


fr © 


ZOO fr *> H <0 


O P- • 


</) X *“"• » w j 


X -1 


ii z u. r: 


*— "t fr 


r ^ «► m* 




— » i < cr 




r h 




<r CV T ~ IV 


u? ^ 


i * e ^ x — 


<UJH* 




r k< k 


1 - frfr CC JC 


x ^ o 


cr << < 


c cr Q •— 


lL : k Lu U. -■ -|L' 


il 3 oei-h 3 : x cr 




frfr CV 


IT 1 IT 


a 


rr sr 


* 


O M 

fr-i 




u 


o o 



o 

ERIC 



100 



r\j nido onMn 

Z Z £ ££!£ 2* £t£: £: <v: tv tv (M nmmn m ro r* ; m ro ro s* 

m mm cnmmmmmmmmmmmmmmmm m m m m mm mmmmmm mm 



it 

3 

CD 

Z 

*-* 

z ~ 
z X 

U w 

r <r 

3T Uf 

a 

3? 5)0 

2Ui 
3 V) 

o 3 

O fr— 

X <t Ui 

u cc 

l/)h <| 

a 

ac < co 

C hUl 
k WH 
3 



CO 

Ui 



o 

UI 

I/) 

3 



3 
O 
X O UI 

uia: 
u <o < 
O 3 

it 

fr- UJ UJ 
UJIH 
it <3 
O 

LlM/) QC 
X u»* 

•— »- 
3 — 
a 

o a ex 



CO CO o 


CO 


UJ Ui X 


O CO 3 


fr- h 


U 2 X 


3 3*** 


3 X X 


OO—i 


C 2 < 


a x x 


3 *~iX 




3 X 


U -ICO 


< 2 


C 3 3 


2 UJ 


< X 


UiUi X 


fr* X 


5X3 


UI o c 


h-t 3 


CO 2 X 


fr- 


< 





< _ 
ou* 

^ UX 
°(^2 

O 

Z UJ£D 
0X3 
C < — 

OD CO -« 

CO < X *-i 
a tu ~ ~ 
ui e- uj cr ~ 

JDIUJO 

o o < a cd 
3 a it 3 -j 
cr 

u K 2 II 

«-* fV' •? «— < 

a a* a 



cm 

o ♦ 

CD 

-J z 



< XI- 
►- X 
O J< 

I- C X 

o 

r x tt 
3 o < 

X it 

M UJ 

XXX 
»-« CV < O «* 

X u <o X 

Ui 

i a 

x 
tv 

fr- II 



O' <r 
a* cr 
o a* 

Ch o 
X* ;> 
^ <* 

• » 

o o 

♦ ♦ 



CM 



« I 

X z 



H II H II 



CM »— f 

or a 
-i 2 



*-« rvx xx 

5 5 <r < u 
ZZX XX 



X 

U 

a 

♦ 

•— * 
— x 
o 3 

• 

II 

ti 

X X 
u 3 

X X 



e> 
2 
U H 



2 CO Q 


U> 


UJ M 


H 


h- < cr 


3 


<r 


O 


a o a 


X 


UI h UJ 




Z Z CD 


it 


llJ M Z 


3 


CD UJ 


X 


3 it 




C U CO 


u 


2 fr- < 


h 


< a. 


►M 


u a 


2 


x ui cr 


•M 


2 •> fr- 


u 


U U > 2 


2 C 


hCh m 


►H O 



3 ~ < 

ouu(r 

(X X < h- 
a w 
i/> i- <i z 
3 0 0 0 
X 2 O 
2 1/5 

U o 3 > 
1-OlEh 

M M 

2 UJ 3 

H K >• M 

U 3 U X 
2 C* CO < 
»r*QCH J 



X 
LJ fr— 
X UJ 
I- X 

cc >» 



fr- o 
<r r 
X h- 

uj <o <o 

2 UJ x 

uj cr x 
cc < »- 

UJ M 



*-« f\J LU 3 UJ UI 

OOI-H-XVO UJ X X 
• • h-h- X X X UJ o 



3 

o 

X 

10 

X 

UJ 



u 

2 



(O 

Ui 

X 

u 

X 

h" 

> 

on 

o 

u 

a 

H 

3 

e 

UJ 

a 



X X 


h »- 


u z O 


u. 


3 3 


h- CO < < 


O O 


X 


X x 


2 3 CO > 


u c 


»-» 




uj a <t 


C X 3 


1- 


* * 


a x 


C fr- < 






ct uj u cr 


X 


c 




3»-i-2 


h- Uj UJ 


z 




O ►-« X flD 


UJ CO CO 


»M 


or er 


2 3 


5' 3 3 


3 


UJ UJ 


u ►-« 




UJ 


a a 


O U <0 2 


X it it 


> 




2 3 2 


0 2 2 


<t 


4 4 


X *- fr- < 


U < <. 


a 



oj ^ x 

OD w w 

-J ^ ^ H 



X 3 U 
2 > 

II UI O 

It 

II II II *-> CV I 
it 



U. M OU *H(V) 



n it n •» 



CV •— (V h- 
X k— h- X 
-J R- K- < 
X X X X 



k cr u 

X <X k- 

c a o 
Z22 



to 

cr 

<r 

a 

z 



nj tv cm 

U U UJ 
3 0 3 
coo 
x x x 



ui 

h* 

3 

O 

a 

CO 

3 

00 

UI 



u 

2 



CO 

Ui 

00 

UJ 

X 

fr- 

ee* 

2 



o • 

o — 

fr- 

cr ii 

JNQ 

<t h*K H O 

h- o 

O Mil II 

B h N > O' 

s£ ~ 

Mil || || 

U w 

O »-*«—• «M 

eh ha c: or 

^ H MM O *4 



UUO 



o oo 



ouu u uuu o u 




101 



*-*cu •HfunsMT^Ncca cvj n * 

Ifi O N CD O' Ort WW^myON CD O C> •— i «~* 

st sf^^^sj-sfrimnirtnintnintninm^NCN©*) *c o >c *> *c v; \C *c %c^> 
men mmnfnrofnnMncoromro 



<Z) 








in 




X 


















U 








3 






















fr- 








Z 




X 


















3 




X 




H 




■-» 


















C 




u 




X 




a 


















q: 




a 




£E 




i- 






















3 




UJ 




< 


















u. 


V) 


tfl 




1- 




X 


















o 


UJ 


2 




























fr- 






Q 




u. 


















►— 


3 


X 




2 




o 


















UJ 


o 


O 




< 






















</> 


a 


X 






























u 




Z 




0"* 


















h 


u 






M 






















<z> 


o 


en 




e 






















tu 




u 




H 




3* 


















cc 


fr- 


fr- 




X 




O 




















UJ 


3 




o 




CC 


















u* 


IT 


a 
















-J 










X 




oc 




e 




o 








O 










fr- 


fr- 






z 




fr- 








c 












lZ) 


u 




H 












X 










>- 


UJ 


o 




a 












o 










CC 


CD 






3 




or 








cz> 


<0+ 












fr- 




_J 




o 














z 




o 


Ui 


UJ 




O 




X 








a 


\D 








UJ 


X 


<z? 




fr- X 












o 






#• 




a 


fr- 






UJ UJ 












u. 






— t 




M 




*- 




\r 




X 


















3 


ee 


IT 




c c 




M 








c 






it 




o 


z 


UJ 




Ui UJ UJ 




X 








fr- 


If 








UJ 


M 


a 




X Oh 




fr— 








< 


»- 


or 


3 


sr 


a: 


fr- 






1— H K. 




<x 








a 


^ a 


c 








< 


u« 




> IZ3 




X 










X 


X 






Ur 


OC 


X 




z a *- 












-j 


• < 






. » 


X 


UJ 


fr- 




•-• UJ > 




X 








-j 


— X 


X 


3 


K 


H- 


Z 






V* 




c 








<3 


»-• 


•~ 




IV 


fr- 


UJ 


ee 




v) in 




a 










*— 4» 


X 


H 




c 


z 




u tr a 




u 










X 


fr— 


— r 




ee 




►- 




fr- a o 














>c o fr* 


< 


a 




X 


IT 


z 


IT S C ►- — 












fr- 


If V X 


5 


c 


c < 


►* 


<Z) 


fr-i 


«~o ►— U) 












3 


-j < 


X 


X 


a *-» 


-J 


c 


< 




a in » 




p-i 








O 


X 


c 




X o 


UJ 


CL 


fr- 


fr- 


u. 




it 








1 


» ^ 




►-4 , — p 


m m 


> 




z 


a: 


u u c — 




w 








fr- 


IV 






< 


u 


o 


2' 


o c •-* ^ 




or 


2 






zc 






•» «K ^ 


X 


o 


o ^ 


fr-^ 




UJ 




3 




*-l 




X X 




X X 


fr— 






XX X Cfl'- • 




X 




#• 




X 


3 4' tn o in o 


o —• 


in m (vj 




a 


X 




U- UJ UJ ►— - V; *-f 




3 ~ 


rx 


* 




X 


► O « O **H 


2 


OJ (V 


(V 


JNUJ 


OlkthX JD 




IZ> 




*** 






MO O O O 


» 


o o 


• » o 


< N. tp 


fr- 


O X 3 X X II O 




z *: 


It 


X 




ii) I! M rH H H 


o » 


I-* r-4 


O H 


t— 


^ r 


o 




3 03 ♦ X 




<««• 




c 




fr— 


IP-4 


X 




X X 


O 


r-3 


UJH 2 QCZ • <T *-> tl 




ar 3 x 




<r 


» * • » 


-4 




#-4 • 


fr- 


t^z 


o > 


< 


• X 




3 uJ 










^ sC «c 


M 


^>40 


>T 


UJ 


ii ii 




X 


II X II ~ 


uj o a 




II 


UJ o 








UJ w ^ 


3' ii 


II II 


cr 


O (V ♦ fHM 


3 


a 3 


h- 




3 UJ ^ ^ *-*■ w 


Hi 


w w 


3 w 








c 


^ C f» fW — 


2 


i/ 


r- 








* 


fr- 




Z fr— fr— 


mk K l/) vlK U 


Ui CZ> 2 ^-r fr- 


M 


fr- z 




3mq: 


UUJUJUJ<lh»UJUJ»-*< C UJ 


fr- fr- 


fr- IT 


IT CL 




ee c. a 


h- U 




•> 


fr- u 




? iT 


fr- fr- 


fr- X X fr- 


x z 


Z <r 


«* 2 


ujct a cr sr ii .z 


Z •-« II 






Z fr- 


•<h v-i >— m cr 


1— I M 


zer xm 


c •- 


*- a 


a *- 


X 


C C C C c h 


c: X 


c 




C 




u 3 a a a 


c c* 


a a 


cccor 


O X X X 


X X fr- S 3 * o X 


OI/> 3C 


QZUm 


H3232 


b. O "* 3 


OIL U. 2 



o 

< 



< 

if 



K C*#- 

if rv rv 

o o 



u o ouooou 



u 



u 




102 



W/? Format ( 1H0* 5 X* bhmatrix M ) 3625 

on 5B I = l* M 3626 

WRITE ( 6* 10?0 ) I 3627 

Write ( 6. 1021 ) ( MMJi, J « 1, N ) 362R 
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irTp c vector containing THf set of routes O0TAINFD FROM The 375 

CURRENT INFINITE Bus ROUTE 376 

The Format OF IRTP is the same AS The FORMAT OF MINRT 377 

UPDATE npass 378 

I9j NpASS S NPASS ♦ 1 379 
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